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The spin-nematic state has proved elusive, due to the difficulty of observing the order parame¬ 
ter in experiment. In this article we develop a theory of spin excitations in a field-induced spin- 
nematic state, and use it to show how a spin-nematic order can be indentffied using inelastic neutron 
scattering. We concentrate on 2-dimensional frustrated ferromagnets, for which a two-sublattice, 
bond-centered spin-nematic state is predicted to exist over a wide range of parameters. First, to 
clarify the nature of spin-excitations, we introduce a soluble spin-1 model, and use this to derive a 
continuum field theory, applicable to any two-sublattice spin-nematic state. We then parameterise 
this field theory, using diagrammatic calculations for a realistic microscopic model of a spin-1/2 
frustrated ferromagnet, and show how it can be used to make predictions for inelastic neutron scat¬ 
tering. As an example, we show quantitative predictions for inelastic scattering of neutrons from 
BaCdV 0 (P 04 ) 2 , a promising candidate to realise a spin-nematic state at an achievable h ~ 4T. We 
show that in this material it is realistic to expect a ghostly Goldstone mode, signalling spin-nematic 
order, to be visible in experiment. 

PACS numbers: 75.10.Jm, 75.40.Gb 


I. INTRODUCTION 

The spin-nematic state is a “hidden order” of spin 
degrees of freedom, involving the ordering of spin- 
quadrupole moments, in the absence of conventional spin- 
dipole order [cf. Fig. 1]. The spin-nematic state was first 
proposed several decades ago 1-3 , and the theoretical pos¬ 
sibility of spin-nematic order is now well-established, es¬ 
pecially for spin-1/2 frustrated magnets in applied mag¬ 
netic field 4-20 . Nevertheless, to date, the spin-nematic 
has generally lived up to its epithet, and remained well- 
hidden from experimental observation. 

The reason why spin-nematic order is difficult to ob¬ 
serve is that its order parameter, a quadrupole moment 
of spin, does not break time-reversal symmetry 3 . This 
means that spin-nematics are invisible to common probes 
of magnetism : they do not lead to magnetic Bragg peaks 
in elastic neutron scattering, asymmetry in muon spin 
resonance (/iSR), or splitting of spectral lines in nuclear 
magnetic resonance (NMR) experiments. In this article, 
we continue the program begun in [21-23] of exploring 
how the symmetries broken by spin-nematic order man¬ 
ifest themselves in its excitations, and how these excita¬ 
tions might be observed in experiment. To this end we 
develop a general theory of inelastic neutron scattering 
from a spin-nematic state in applied magnetic field. 

The scenario we explore, summarised in Fig. 1 and 
Fig. 2, is applicable to a wide range of materials. When 
a frustrated magnet is polarised by applied magnetic 
field, interactions between magnons can lead to the for¬ 
mation of a two-magnon bound state. At lower val¬ 
ues of magnetic field, this bound-state can condense, 
leading to spin-nematic order 4,5,24,26 [Fig. 1]. Since the 
spin-nematic order breaks spin-rotation invariance in the 



FIG. 1: (Color online). Two-sublattice, bond-centered spin- 
nematic state, of the type found in spin-1/2 frustrated ferromag¬ 
nets in applied magnetic field. Spin fluctuations show quadrupolar 
character, visible in the probability distribution for spin fluctua¬ 
tions on the bonds of the lattice, here represented by a blue sur¬ 
face. Spin-nematic order is known to occur close to saturation in 
the spin-1/2 J 1 -J 2 model on a square-lattice, [Eq. 1], for 

ferromagnetic J\ and antiferromagnetic J 2 [5,9,24]. This model is 
believed to describe a number of quasi-two dimensional magnets, 
including BaCdV0(P0 4 )2 [25]. 

plane perpendicular to the magnetic field, it must pos¬ 
sess a Goldstone mode. This has observable consequences 
— the two-magnon bound state, invisible to neutrons in 
the polarised phase, transforms upon condensation into 
a ghostly, linearly-dispersing Goldstone mode, which can 
be resolved in inelastic neutron scattering [Fig. 2]. 

For concreteness, in this article we concentrate 
on spin-1/2 frustrated ferromagnets on a square lat¬ 
tice, taking our motivation from materials such as 
Pb 2 V0(P0 4 )2 27 “ 31 and SrZnVO(P0 4 )2 29 ’ 30 ’ 32 . We pay 
particular attention to the quasi-two dimensional frus- 
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FIG. 2: (Color online). Predictions for inelastic neutron scattering from a powder sample of a material exhibiting bond-centered spin- 
nematic order in applied magnetic field, of the type shown in Fig. 1. (i) Inelastic scattering at small |q|, for magnetic field approaching 
the saturation value h sa t- The existence of a spin-nematic state is heralded by a ghostly, linearly-dispersing Goldstone-mode at low energy, 
(ii) Inelastic scattering at fixed |q| = 0.18A -1 , showing the distribution of spectral weight as a function of frequency. The Goldstone mode 
has at its maximum about 3% of the intensity of the spin-wave mode at higher energy. Predictions for the dynamical structure factor 
X°^(q, w) were obtained using the methods described in Section VI of this article, for parameters relevant to BaCdVOfPCqE [25], and 
powder-averaged. Equivalent results for a single-crystal sample are shown in Fig. 16. 


trated magnet BaCdV 0 (P 04)2 25,29 , which has a satura¬ 
tion field of only h sat ~4T, and so is easily accessible to 
experiment. However, with small modifications, the same 
methods and conclusions can be generalised to other sys¬ 
tems, such as coupled spin-half chains 16-18,20 

The primary microscopic model we consider is the 
spin-1/2 “J 1 -J 2 ” Heisenberg model, 

=Ji E s '- s ./ + J 2 E s - s i + h E ^ fi) 

di> l (ij) 2 i 

where (ij) i counts first-neighbour bonds, and (ij) 2 
second-neighbour bonds of a square lattice [see Fig. 1]. 
Magnetic field h defines the S z axis for spin. This model 
is believed to describe several distinct families of quasi- 
2D materials, including the vandates Pb 2 V 0 (P 04 ) 2 27-31 , 

SrZnV0(P0 4 )2 29, 30,32 and BaCdV0(P0 4 ) 2 25,29 . 

The “J 1 -J 2 ” model [Eq. (1)] can be shown 

to support a spin-nematic ground state for all ferro¬ 
magnetic J\ < 0 and antiferromagnetic J 2 > 0.4081 Ji| 
[24,26]. Spin-nematic order is formed through the con¬ 
densation of bound pairs of magnons out of the satu¬ 
rated state at h = h sat [5], and is generally believed to 
be stable for a small range of fields approaching satura¬ 
tion, i.e. for h < h sat . However for a range of parameters 
0.408 < J 2 /I J\\ < 0.7, the zero-field (h = 0) ground state 
of [Eq. (1)] is also a spin-nematic state 5,9,33-35 . 

Id and quasi-Id frustrated ferromagnets have also 
been extensively studied in the search for the spin- 
nematic state. Theoretically it has been shown that 
J 1 -J 2 , spin-1/2 chains in applied magnetic field, with 
Ji < 0 and J 2 > 0, demonstrate dominant quadrupolar 


correlations for a wide parameter range 4,7,8,10-15, 36-38 # 
In the presence of small interchain coupling this can 
lead to a long-range-ordered spin-nematic state at low 
temperature 16 18,20 . The spin-nematic state is stabilised 
by magnetic field, and is most pronounced close to 
the saturation held. There have been a number of 
calculations of dynamical properties of such a spin- 
nematic state, with a view to providing experimen¬ 
tal predictions 16-18,20,39,40 . The material LiCuV 04 is 
thought to approximately realise this model, and may 
show spin-nematic order close to saturation 16,41,42 . How¬ 
ever, high-held NMR measurements have not yet de¬ 
tected evidence for such a state, and have shown that, if 
it does exist, it is limited to a very narrow held range 22,42 . 


The spin-nematic state that appears in spin-1/2 mod¬ 
els such as [Eq. 1] is known as a bond-nematic. 

While an individual spin-1/2 cannot have a quadrupole 
moment, a pair of neighbouring spin-l/2’s can form a 
triplet, and thus develop a quadrupole moment living 
on the bond 3 . If conventional dipole magnetism is sup¬ 
pressed, for example due to high frustration, quadrupo¬ 
lar or higher-order multipolar correlations can be re- 
vealed. In j~Lj~_j [Eq. 1] the triplets organise them¬ 
selves into a bond-centered antiferroquadrupolar (AFQ) 
order 5,9,33-35 . The way in which this occurs is most easily 
understood at high values of magnetic held, in the satu¬ 
rated state, where triplets are preformed. Here magnons 
form bound states, which condense to give spin-nematic 
order as the magnetic held is lowered 5 . In such a state 
the spin-dipole moment (S) = 0, while the rank-2, sym- 
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FIG. 3: (Color online). Detail of a bond-centered, two-sublattice 
spin-nematic state in applied magnetic field. Green spheres repre¬ 
sent spin-1/2 degrees of freedom at the vertices of a square lattice. 
The bond-centred, nematic order parameter is represented by red 
cylinders, and the associated distribution of spin fluctuations by 
a blue surface. A second, bond-centred lattice is introduced and 
shown in yellow. In this article we consider both the site-centered 

S=1 /2 

spin-1/2 Heisenberg model j 2 [Eq. 1], and a bond-centered 
spin-1 bilinear-biquadratic (bbq) model [Eq. (4)]. Both ex¬ 

hibit the same form of spin-nematic order, in applied magnetic 
field. 

metric, traceless tensor, 

QTf = s?sf + S?S? - jj<5 a/3 S i ■ Sj, ( 2 ) 

with a,/3 = x, y, z, has entries with non-zero expectation 
value 3 . 

At h = 0, one possible approach to understanding the 
excitation spectrum of the spin-nematic state in %ji_j 2 
[Eq. 1] is to work directly with the microscopic model, 
construct a lattice gauge theory, and then solve this 
within a large-N expansion scheme 33-35 . This approach 
has the advantage that, in principle, it can access all 
of the different excitations of the spin-nematic state. 
However, the solution of the lattice gauge theory is ex¬ 
tremely involved, which complicates the interpretation of 
experiments 35 . 

A second possibility — explained in detail in Refs. [21, 
23] — is to construct a continuum theory for the long- 
wavelength excitations of 2-sublattice AFQ order. This 
approach has the advantage of bringing the universal 
properties of the spin-nematic state to the fore, and of 
making clear predictions for inelastic neutron scattering. 
However, being grounded in the symmetry of the order 
parameter, it cannot hope to describe the microscopic de¬ 
tails of the underlying spin-1/2 model at high energies. 

In this article we combine the continuum theory ap¬ 
proach with a microscopic study of J-Ln-n [Eq. 1]. We 
use diagrammatic calculations to determine the magnetic 
dispersion spectrum of 1-magnon and 2-magnon excita¬ 
tions at and just below h sat . This allows the continuum 
theory to be parametrised using J\ and J 2 , and therefore 
quantitative predictions to be made for inelastic neutron 
scattering experiments. 


Inelastic neutron scattering measures the dynamical 
spin correlation function, which is defined as, 

Amy a/ 3 (q, uj) = 

POO 

/ dte™(8SZ(t)SSt^ 0)}}, (3) 

Jo 

where we have set g/is = 1. Thus the task of this article 
is to calculate this quantity, and we will do this for two 
complementary models of field-induced, spin-nematic or¬ 
der on the square lattice: Un-n [Eq. 1] and a spin-1 
model that is a generalisation of the bilinear-biquadratic 
(BBQ) Hamiltonian 1,2 . We now briefly review the route 
taken. 

In Section II we introduce a spin-1, BBQ model with 
a partially polarised, 2-sublattice, spin-nematic ground 
state. As shown in Fig. 3, one way to view this is as an 
effective model describing spin-1 degrees freedom living 
at the bond centres of a spin-1/2 square lattice. The lat¬ 
tice of bond centres also forms a square lattice, with a 
reduced lattice constant b = a/y/2. A major advantage 
of spin-1 models is that the excitation spectrum can be 
calculated within flavour-wave theory 43-47 . This allows 
predictions for Srnx a/3 (q, uj) [Eq. 3] to be made. While 
the primary motivation for considering this model is as 
a first step towards making predictions for spin-1/2 sys¬ 
tems, spin-1 is also interesting in its own right 46,47 , and if 
real materials can be synthesised with large biquadratic 
interaction, the results presented here would be relevant. 

In Section III we start from the spin-1 BBQ model 
studied in Section II and use it to derive a continuum 
field theory. We first demonstrate that, at long wave¬ 
length, this exactly reproduces the flavour-wave results 
of Section II. The power of the field theory is that it is 
a theory of the order parameter symmetry, and therefore 
describes the universal features of a partially polarised, 2- 
sublattice spin-nematic state. Thus we recast the theory 
in terms of a minimal set of hydrodynamic parameters. 
This renders the theory free of any particular microscopic 
model, and one can in principle parametrise it from any 
microscopic model with a partially polarised, 2-sublattice 
spin-nematic ground state or directly from experiment. 

In Section IV we make a mapping between the effective 
spin-1 degrees of freedom on the bonds and spin-1/2 de¬ 
grees of freedom on the sites. This will allow predictions 
to be made for spin-1/2 frustrated ferromagnets, based 
on the theory developed in Section II and Section III. 

In Section V we consider 'H Ji Tj 2 [Eq. 1] from a mi¬ 
croscopic perspective close to the saturation field. We 
consider the condensation of magnons out of the fully- 
saturated state. For a sizeable parameter range, conden¬ 
sation of bound-magnon pairs occurs at a higher mag¬ 
netic field than the condensation of single magnons, and 
therefore a spin-nematic state is formed. Diagrammatic 
calculations allow the critical field to be calculated, as 
well as the velocities and gaps of the excitation modes 
in the spin-nematic state. These can then be used to 
parametrise the continuum theory. 
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In Section VI we make predictions for inelastic neu¬ 
tron scattering experiments for materials described by 
[Eq. 1]. As a worked example we consider the 
material BaCdV 0 (P 04 ) 2 , which is expected to have a 
spin-nematic ground state close to saturation, and show 
quantitative experimental predictions. 

Finally in Section VII we conclude by showing 
that the detection of a ghostly Goldstone mode in 
BaCdV 0 (P 04)2 - a characteristic signature of spin- 
nematic order - is experimentally feasible using current 
instruments. 

In Appendix A we derive a non-linear-sigma-model 
field theory for the 2 sublattice AFQ state at h = 0. This 
is complementary to the continuum theory presented in 
Section III, which considerably simplifies at h = 0 if a 
set of high-energy modes are eliminated by a Gaussian 
integral. This theory can then be compared to previous 
work considering the spin-nematic state of [Eq. 1] 

at h = 0 33-35 . 


II. A SPIN-1 MODEL FOR THE 2-SUBLATTICE 
SPIN NEMATIC IN APPLIED MAGNETIC 
FIELD 


Here we construct and solve a spin-1 bilinear- 
biquadratic model in applied magnetic field on the square 
lattice that supports the same type of spin-nematic state 
as is found in [Eq. 1]. While this may be rel¬ 

evant to spin-1 systems, the primary motivation is as 
an effective model of bond degrees of freedom in a spin- 
1/2 bond-nematic state. We make the assumption that 
the 2-sublattice, spin-nematic state reduces to triplets on 
nearest-neighbour bonds and, at low energy, singlet de¬ 
grees of freedom can be ignored. The exchange parame¬ 
ters in the spin-1 model are chosen such that a partially- 
polarised, 2-sublattice AFQ ground state is realised for 
the full magnetic field range 0 < h < h sat . 

We use flavour-wave theory 43-47 to determine the evo¬ 
lution of the magnetic dispersion and the imaginary part 
of the dynamic spin susceptibility as the magnetic field is 
varied. Here we present the results of linear flavour-wave 
theory (see Fig. 6). At h = 0 we have checked that the 
2-sublattice AFQ ground state remains stable when in¬ 
teractions are included, but the details will be presented 
elsewhere 48 . 


A. Definition of the model 

The model we introduce is a straightforward general¬ 
isation of the bilinear-biquadratic model introduced in 
Ref. [1,2]. Effective spin-1 degrees of freedom are placed 
at the bond centres of a square lattice with lattice con¬ 
stant a. The bond centres of the square lattice also form 
a square lattice, with lattice constant b = a/yj 2 (see 
Fig. 3). 


The microscopic model is given by, 

n^[s] = E j n [(Si.s,-) 2 + s^] -h^si 

(ij )i i 

- E [ J 12 [(Si-S,) 2 + S,.S,] + J 22 (s i .s i ) 2 ] , (4) 

(ij) 2 


where (ij) i counts the first- and (ij )2 the second- 
neighbour bonds on the bond-centred lattice [see Fig. 3] 
and S «= (S x ,S y ,S z ) is the usual spin-1 operator. We 
consider all interactions J to be positive and h is the 
applied magnetic field. 

It is useful to rewrite ^^[S] [Eq. 4] as, 


ni^[s, Q] = E w + s * ,s i) ~ h Y, s i 


-E 

hi) 2 L 


di) 1 i 

~ (Qi-Qj + Si.Sj) + j -f (Qi.Qj - Sj.Sj) 


( 5 ) 


where, 


Q = 


Q x2 - y2 \ 

Q xy 
Q yz 

Q xz ) 


(S x ) 2 - (S y ) 2 \ 

-L[2(S*) 2 - (S*) 2 - (S y ) 2 } 

s x s y + s y s x 
s y s z + s z s y 
s x s z + s z s x ) 


(6) 


describes spin-quadrupole operators and a constant term 
has been dropped. 

This model has a hidden SU(3) symmetry for J 22 = 0 
and h = 0, a fact which is more easily understood if 
Hf^[S] [Eq. 4] is expressed in terms of a director vector, 
d. This follows from noting that the wavefunction of a 
spin-1 on a site j can be written as 21,44,49-51 , 


\d j ) = d x \x)+d y \y) + d z \z), (7) 


where, 


i*> = 


|i) + |i) 
72 


17 


-70), 


(8) 


are linear superpositions of the usual spin-1 basis states. 
The director vector is defined by d j = (d*,d 7 ,d*) and is 
normalised by requiring dj.dj = 1. Assuming the total 
wavefunction can be site-factorised allows ^bbq 1 [S] [Eq. 4] 
to be re-expressed as, 


^bbtld] = Jn E I- J 12 E l d -4f 

(ij) i (ij) 2 

- J 22 E l d i- d i! 2 - i/l E ~ d ^ X i) • 

(ij) 2 i 

The first and second terms of ^bbq 1 [d] are SU(3) invari¬ 
ant, and favour general two-sublattice states. However, 
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h=h sat h=0.75h sat h=0.5h sat h=0.25h sat h =0 


FIG. 4: (Color online). Evoluiton of the 2-sublattice, AFQ spin-nematic state in applied magnetic field. The distributions of spin 
fluctuations, shown by blue surfaces, are calculated within the mean-field ground state described by Eq. 15. Magnetic field is applied 
perpendicular to the 2-d plane, and polarises the magnetic moments. Full polarisation is achieved at h = h sa t- The nematic order parameter 
is shown by red cylinders, and disappears at h — h sat . In the case of a spin-1 system, the polarised moments and the quadrupolar order 
parameter exist at the vertices of the yellow lattice. In the case of a spin-1/2 system, magnetic moments are associated with the vertices 
of the green lattice (shown by green spheres). The nematic order parameter lives on the bonds, and one can define partially polarised 
spin-dipoles on the bonds by summing contributions from neighbouring sites (see Section IV). Magnetic field values correspond to those 
in Fig. 6 and Fig. 7. 


the interaction J 22 breaks the symmetry of the model 
down to SU(2), and, enforces two-sublattice AFQ order 
of the type shown in Fig. 3. The magnetic field further 
reduces the symmetry to U(l) and, above a critical field 
ft- sa t, favours a saturated paramagnetic state. 

An alternative rewriting of T^bq 1 [S] [Eq. 4] is in terms 
of a vector e, which parametrises the usual spin-1 basis 
states. The wavefunction on a site j is written as, 

l e j) = e j|l) + e°|0) + ej|l), (10) 

and Qj = (ej, e°, ej). The wavefunction is normalised by 
requiring ej.ej = 1 . Again assuming a site factorisation 
of the wavefunction, ^bbqH^] [Eq. 4] can be rewritten as, 

= Ji 1 J 2 l e *‘®il 2 “ J i2 E 

W) 1 (ij) 2 

- J22 E \ e l e ] - e °i e °j + e i e )\ 2 -^E(ldl 2 - Idl 2 ) • 

(ij) 2 i 

(ii) 

At large values of h is it clear that all sites will have 
e == (i, 0, 0), corresponding to a saturated paramagnet. 

The relationship between the e vectors and the opera¬ 
tors S and Q is given by, 


B. Mean-field ground state 

The mean-field ground state of ^bbqM^] [Eq. 4] can 
be derived by varying the e vectors in ^bbq 1 [ e ] [Eq* H]- 
Since the motivation for considering the model is to study 
the partially polarised spin-nematic state, we do not 
present the entire phase diagram, but instead determine 
a parameter range in which this state is stable for the 
full field range 0 < h < h sat . This is the case for Jn = 1, 
J 12 > 0, J 22 > 1, and from now on we concentrate on 
this region in parameter space. 

In the saturated paramagnet the mean-field wavefunc¬ 
tion at every site is described by, 

e™ f = i (1,0,0). (13) 

Below the saturation magnetic field, 

^sat = 4(Jn + J 22 ), (14) 

the e vectors cant, forming a 2-sublattice state. Labelling 
these two sublattices A and B, the mean field ground 
state can be described by, 

ef = R(0 h ) • ef, eg f = R(-0 h ) • ef, (15) 

with, 

( cos Oh 0 sin 0ft, \ 

0 10. (16) 
— sin Oh 0 cos Oh J 


(S x ) \ 


<s y > 


<5 Z > 


<e x2 - y2 > 


<e 3z2 - r2 > 


(Q xy ) 


(Q yz ) 


(Q xz ) J 



( — ~^( ele ° + ^e 1 + e°e 1 + e 1 e°) ^ 


V2 


(e 1 e° — e°e 1 + e°e 1 — e 1 e°) 
„112 \A\2 


e 1 e 1 + e 1 e 1 


7sd el l 2 + |e 


e 

1|2 


V2 

V 75(- elg ° 


2|e°| 2 ) 

— 6^6^) 

(e 1 e° - e^e 1 - e 0 e T + eh 0 ) 

e^e 1 + e°e l + e^e 0 ) J 
( 12 ) 


The canting angle, 6^^, is given by, 

cos 2 0 h = A ( 17 ) 

hsat 

For h = h sat the canting angle is Oh =0, while for h = 0 it 
is Oh = 7r/4. The field evolution of the mean-field ground 
state described by Eq. 15 is shown in Fig. 4. 


C. Linear flavour-wave theory in the saturated 
paramagnet 


where (S x ) = (e| *S rx |e). 


Above a critical magnetic field, h sat [Eq. 14], the spins 
align, forming a saturated paramagnet. The excitation 
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FIG. 5: (Color online). Relationship between the site-centred lattice of _^j 2 [Eq. 1] and the bond-centred lattice of [S] [Eq. 4] 
and associated Brillouin zones (bz). (i) The site-centred lattice, with lattice constant a, is shown in black and the bond-centred lattice, 
with lattice constant b = a/V2, in red. A coordinate system (x x ,Xy) is associated with the sc lattice and (r x , r y ) with the bond-centred 
lattice, (ii) The site-centred bz is shown in black and the be bz in red. The associated coordinates are (p x ,Py) and (g x , Qy)- High symmetry 
points in the site-centred bz are pr = (0,0), px = (7r/a, 0) and pm = (7r/a, 7r/a) and in the bond-centred bz qr = (0,0), qx = (7t/6, 0) 
and qM = (7t/6, 7t/6). 


spectrum for h > h sat can be calculated using linearised 
flavour-wave theory 43 47 . From this the imaginary part 
of the dynamical spin susceptibility can be determined. 

The spin operators are written in terms of a pair of 
boson creation and annihilation operators, labelled a and 
6, according to, 

S j = 1 - 2 a] aj - bjbj 

S+ = -V2 a)aj - b]bj b j + b)a^j 

S~ = ~s/2 (b)\j 1 - a]aj - b^bj + , (18) 

where S± = S x =b iS y . The operator creates an exci¬ 
tation with A Sj = 1, while Oj creates an excitation with 
A Sj = 2. It follows from Eq. 6 and Eq. 18 that, 


Q x y = a) \/l — of a — b^b + \J\ — at a — b^b a 

Q xy = i ^ y/l — eda — b^b — \/l — at a — b^b a ^ 
Q yz = —j= (^\/1 — eda — tfb b 

—b^y/l — cd a — b^b + a^b — b^a^ 

Q xz = -~^= (V 1 - at a - 6+6 b 

+id\/ 1 — a^a — fot b — adfr — 


(19) 


The dispersion of magnetic excitations can be calcu¬ 
lated within linear flavour-wave theory. This involves 


rewriting Hj^bq^S, Q] 5] in terms of the boson oper¬ 
ators defined in Eq. 18 and Eq. 19, and retaining only 
terms up to bilinear order. The resulting dispersion re¬ 
lation has two branches, 

= - 4Jn(l - lF) + 4 Ji 2 (1 - lF) 

-4J 22 (1+ 7 F) + 2/i 

w k,/i = _ 4Jn(l - 7 F) + 4 Ji 2 (1 - 7 F) + h (20) 
where, 

7 ^ = -(cos k x + cos fc y ), 7 ^ = cosk x cosk y . ( 21 ) 

The imaginary part of the dynamical spin suscepti- 
bilty, Smx a ^(q, cj) [Eq. 3], can be calculated within the 
flavour-wave approach, yielding, 


Smx xx (q, w) = Qmx m (q,oj) = ^S(oj - w,) 
^smx zz { q, cj)=0. (22) 

The dominant feature is a band of 1-magnon excitations 
in the perpendicular channel, as shown in Fig. 6 (a) and 
(b). This is gapped at all wavevectors, and has uniform 
spectral weight. There is also a band of 2-magnon ex¬ 
citations in the longitudinal channel with zero spectral 
weight. For h > h sat this is gapped at all k, but exactly 
at h — h sat the gap closes at k = 0. This signifies the 
onset of spin-nematic long-range order. The spin-dipole 
remains parallel to the applied magnetic field, and a di¬ 
rector order parameter appears in the plane perpendicu¬ 
lar to h [see Fig. 3]. 
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D. Linear flavour-wave theory in the partially 
polarised spin nematic 


where 7 ^ and 7 ^ are defined in Eq. 21. The second 
branch has, 


The dispersion of magnetic excitations in the par¬ 
tially polarised spin-nematic state, which occurs at 
0 < h < ft, sat , is now calculated using linear flavour-wave 
theory. This is done by considering small fluctuations 
around the mean field ground state given in Eq. 15. 

The spin operators can be rewritten as, 


S x = — —(cos Oh — e lku ' Vi sin Oh) ^71 — ata — b^b b 
+&t\/1 — at a — &t b'j + (cos Oh + e zkM ' ri sin Oh) (o)b + b^a) 
S\ = ~^= (cos Oh + e 2kM ‘ r * sin 6 ^) ^\/1 — ata — b^b b 

—b^\/ 1 — ata — tfb'j + (cos Oh — e zkM ' ri sin 6 ^) ( b^a — o)b) 
= cos 2 Oh (l — 2 a^a — b^b) 

+ e zkM ' ri sin 20 h ^7 1 — ata — tfb a + at 7l — ata — frtfr^ 

(23) 


while the quadrupole operators are given by, 


2 2 
<% _y 

= - 

. e «k M -ri 

sin 2 Oh (l - 

• 2ata — fttfr) 



+ cos: 

2 9 h ( 

7i- 

a) a — b^b a 

+ at ^1 — a 

ta- 


Q?~ r ' 


- 3 tfb) 




II 

x 

Q 1 

i 

Vi - 

at a — fotft — 

7l — at a — 


) 

II 1 

N 

Q 1 

i 

71 

(cos 0 t 

h — e ik M-ri gin ^) (/[Z 

ate 

t — fttft b 

-6V1 

— ob 

a — b^bj + (cos Oh 

+ e lkM ' ri sin 

Oh) 

(at b — b^a 

Qf = 

1 

71 

— (cos 

6» fe +e ikMri 

sin^) (7l 

— a 

ta — fotfo b 

+bW 1 

-at 

a — frtfrj -f (cos Oh 

_ gtkM-r; g j n 

Oh) 

(at b + ftt a 


(24) 


where kM = (7r, 7 r) is the 2-sublattice ordering vector [see 
Fig. 5] and e zkM ' ri = ±1 depending on whether i is within 
the A or B sublattice. 

Performing linear flavour-wave theory results in a mag¬ 
netic dispersion with two branches. The first has the 
excitation spectrum, 

,h = y/ A i,h B l,h> ( 2 ^) 

with, 

Ak,h = 4Jn sin 2 29 h - 4Jn(l - j^) cos 2 28 h 

+ 4Ji 2 (1 - 7^ 2) ) + 4J 2 2 sin 2 2 9 h - 4J 22 (1 + j[ 2) ) cos 2 20 h 
+ 2 h cos 20h 

Bk,h = -4Ju7k 1) sin2 2 ^h + 4J 2 27k 2) sin2 ZOh, (26) 


w, 


k ,h 


— 7, (Ck ,h ~ Ck+k M ,/») 


+ ^ v /(C k , ?l + Ck + k M) / l ) 2 -4D 2 j „ (27) 


where, 

Ck,h = - 4 Jii cos 2 26^ + 4 Jn 7 k X) cos 

+ 4Ji 2 (1 - 7 ^ 2 - ) ) + 4J 22 sin 2 2 #^ + h cos 29k 
D Kh = - AJ 22l (2) sm28 h . (28) 


It follows that the imaginary part of the dynamic spin 
susceptibility is given by, 

7s : mx ± (q, uS) = Amx xx (q, cj) + 7s7nx yy (q, uS) = 

7T (tiq >h COS 8 h - v\ h sin 9 h ) 2 S(u - w,) 

+ * ( u q ,h s[ n9 h - V^ h cos 9 h ) 2 5(u - w q +q M ) 
9mx zz (q, lo) = i r sin 2 28^ 

( U q+qM ,h + W q+q M - W q+q M )- ( 29 ) 

where, 


«7 2 

^ u ii,h v t,h 


^4k,/i 

b 

/^k,/i ^k,/i 


Ak,/z 

2 \ 

/^k,/i “ ^k,/i 


~B\<L,h 


\f A i,h ^k ,h 


(30) 


«J 2 = 
«,) 2 = 
^ u k,h v t:,h = 


_Ck,fr + Ck+k M ,/i_ 

^k+k M ,h) 2 ~ 

_Ck,/i + Ck+k M ,fr_ 

2 ^/ (Ck,/i + Ck+k M ,h ) 2 — 4-Dj^ 



1 

2 


_—2T>k,/i_ 

^ (Ck,/i + Ck+k M ,/i ) 2 — 


(31) 


The imaginary part of the dynamical spin susceptibility 
[Eq. 29] is plotted in Fig. 6 . An animated version of Fig. 6 
with continuously varying magnetic field is provided as 
supplementary material 52 . 

For h > h sa t [Fig- 6 (a)] the only visible feature is a uni¬ 
formly bright band of A S z = 1 excitations in the trans¬ 
verse channel. The only other excitations of the satu¬ 
rated paramagnet have A S z = 2 and therefore do not 
contribute to the spin susceptibility. 

At h = h sa t [Fig. 6 (b)] the gap to these A S z = 2 exci¬ 
tations closes, signalling the onset of spin-nematic order. 






































8 



(ii) 3mx zz (q, cd 









q 


q 


FIG. 6: (Color online). Flavour-wave predictions for the imaginary part of the dynamic spin susceptibility of a spin-1, partially-polarised, 
2-sublattice, spin-nematic state in applied magnetic field. These follow from considering ^bbq 1 ^] Pq- 4] with Jn = 1, J ±2 =0.1 and 
J 22 = 2. The external magnetic field is gradually reduced from a) h = 1.5h sa t to f) h = 0. (i) The transverse susceptibility Smx 1 (q, co) 
[Eq. 29]. Dashed red lines show a; 5 h at all h and ^q_bq M h ^ or ^ — ^ sat [ see ^q. 20, Eq. 27]. The associated intensity is shown by the 
colour scale inset in panel (ii)a. (ii) The longitudinal susceptibility Srnx zz (q, uj) [Eq. 29]. Dashed red lines show ^q+ qM h [ see Eq. 20, 
Eq. 25]. The associated intensity is shown by the colour scale inset in panel (ii)a. All predictions have been convoluted with a gaussian 
to mimic experimental resolution. The circuit T-X-M-T in the bond-centred Brillouin zone is shown in Fig. 5. An animated version of this 
figure is shown in the supplemental material 52 . 
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For 0 < h < h sat it is no longer possible to assign an 
integer value of A S z to a particular excitation. In con¬ 
sequence spin and quadrupolar fluctuations are mixed, 
except for at special points in the Brillouin zone. Of par¬ 
ticular interest is the nature of the low-energy modes, 
and we consider each of these separately. 

The spin-nematic phase possesses a gapless Goldstone 
mode, which can be seen at q = (0, 0) in the longitudinal 
spin susceptibility, Smx zz (q, ui) [Eq. 29, Fig. 6(ii)(b-f)]. 
This is due to spontaneous breaking of the U(l) symme¬ 
try of H^bq 1 [S] [Eq. 4] in field by the spin-nematic state. 
The Goldstone mode is associated with rotations of the 
quadrupoles in the plane perpendicular to the field, and 
shows up at q = (n, i r) in the h [Eq. 25] branch of 
the dispersion. Close to q = ( 7 r, 7 r) the predominantly 
quadrupolar fluctuation induces a small spin fluctuation. 
These spin fluctuations are parallel to the z direction and 
in phase between the two sublattices, and hence are seen 
at q = (0,0) in the spin susceptibility. Their magnitude 
goes to zero approaching q = (0,0), which can be seen 
from the fact that the spectral weight in the spin suscepti¬ 
bility disappears linearly with |q| approaching q = (0, 0). 

The longitudinal spin susceptibility, ^mx zz (q, cj) 
[Eq. 29, Fig. 6(ii) (b-f)], also has a gapped mode at 
q = ( 77 , 7 r). This is brightest at h = 0 and gradually fades 
as the field is increased towards h = h sat . The mode can 
be thought of as a dynamical spin-density wave, in which 
spin fluctuations occur parallel to the field direction, but 
in antiphase between the two sublattices. Intriguingly, 
a number of quasi 1-dimensional materials that are good 
candidates for realising the spin-nematic state show spin- 
density wave order over large field ranges 41 . It would be 
interesting to see if at different values of the interaction 
parameters this mode condensed. 

The transverse spin susceptibility, ^mx ± ( q, cj) 
[Eq. 29, Fig. 6(i) (b-f)], has a single excitation branch at 
h = h sa t- As the field is reduced below h = h sat it splits 
into two, corresponding to h and <^q +qM ^ [Eq. 27]. At 
a given q there is a bright feature at c 0 = h: which 
is predominantly due to a transverse fluctuation of the 
partially polarised moment. There is also a weaker fea¬ 
ture at c jJ = co q+ qM ft, which is due to the quadrupolar 
part of the fluctation at q = q + qM inducing a spin 
fluctuation at q. Quadrupolar fluctuations at q also in¬ 
duce spin fluctuations at q, but, unless h = 0, these 
are hidden by the transverse fluctuation of the polarised 
moment. At h = 0 there is no polarised moment and 
all spin fluctuations are induced from quadrupole rota¬ 
tions. At h = 0 there are Goldstone mode excitations 
at q = (7r,7r) and q = (0,0) in the transverse chan¬ 
nel, making a total of 3 Goldstone modes in the system. 
The spin-nematic state spontaneously breaks the SU(2) 
symmetry of ^^[S] [Eq- 4], and the Goldstone modes 
correspond to 3-dimensional rotations of the quadrupole 
order parameter [see Appendix A]. 


III. CONTINUUM APPROACH TO THE 
SPIN-1, 2-SUBLATTICE SPIN NEMATIC IN 
APPLIED MAGNETIC FIELD 

Here we develop a continuum approach to understand¬ 
ing the fluctuations of the spin-nematic state. We derive 
a continuum model directly from the microscopic spin- 
1 model [S] [Eq. 4] . After linearising the quantum 
fields, we show that this exactly reproduces the flavour- 
wave results presented in Section II at long wavelength. 
The advantage of the continuum approach is that it is not 
tied to a microscopic model, but is instead a theory of the 
order parameter symmetry. Thus it can be parametrised 
from any microscopic model that supports a 2-sublattice 
spin-nematic state. In Section VI we will parametrise it 
from [Eq. 1], 


A. Continuum theory 

A continuum field theory can be derived directly from 
^bbq 1 [ e ] [Eq- 11]. The spin states on a site are repre¬ 
sented using spin-coherent states, and the overcomplete- 
ness of the spin-coherent state basis results in a geomet¬ 
rical phase term in the action. Assuming that the system 
at least locally realises a 2-sublattice, partially-polarised, 
spin-nematic state, ^bbq 1 [ e ] [Eq. 11] can be expanded in 
terms of derivatives of a set of continuum fields, which are 
defined at the centre of square plaquettes. These fields 
are collected in the complex vectors eA(r, r) and ee(r, r), 
which are the continuum versions of ej [see Eq. 10] on 
the A and B sublattices. The derivation follows in spirit 
Ref. [21], in which a continuum theory is formulated for 
a 3-sublattice spin-nematic state on the triangular lattice 
in the absence of magnetic field. 

The action for the 2-sublattice, partially-polarised spin 
nematic is given by, 

<?2Sl =2 J dtd 2r £2SL, (32) 

where we set the lattice constant b = 1 and, 

£ 2S l = 4sl - (33) 

The kinetic term is given by, 

4sl = ie A .<9 t e A + *e B .dte B , (34) 
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and the Hamiltonian term can be written as, 

£ 2 SL = 4Jn|e A • e B | 2 - 2 J 2 2 | 2 e^ 1) e^ 1) - (e^ 0) ) 2 | 2 
- 2J 22 |2e^ 1) e^ 1) - (e^ 0) ) 2 | 2 - h (|e^| 2 - |e^| 2 ) 


- Jn [(eA • e B )( 9 Ae A • 9 A e B ) + (e B • e A )( 9 Ae B • c^a)] 
+ 2 J 12 [( 9 a e A • 9 a e A ) + ( 9 A e B • 9 Ae - B )] 


+ J 22 [ 24 1 ) 4 1) - (e^ ) ) 2 ][2aAe^ ) a A e^ 1) 
+ J 22 [24 1) 4 1) “ (eA 0> ) 2 ] [ 2 a A e^ L) 0 A e^ 1) 

+ </22[2eg 1) eg 1) - (egVlP^AeB^Aeg^ 

+j 22 [24 1 ) 4 1) - (4 0) ) 2 ][2«a4 1) ^a4 1) 


(3a4 0) ) 2 ] 

(3Aei 0) ) 2 ] 

(Sa4 0) ) 2 ]} 


It is understood that the constraints |e A | 2 = 1 and 
|e B | 2 = 1 have to be enforced. 


When working at h = 0, it is natural to divide the fluc¬ 
tuations into two sets, with one set describing Goldstone 
modes and other low-energy fluctuations, while a conju¬ 
gate set describes high-energy fluctuations. The action 
considerably simplifies if the high-energy fluctuations are 
integrated out 21 [see Appendix A]. However, for fields 
close to h = h sat , it is no longer possible to partition the 
fluctuations in this way. Conjugate pairs of fluctuations 
become degenerate at h = h sat , and therefore there is no 
low-energy fluctuation to integrate out. 


I 

B. Linearising the continuum theory 


In order to calculate the dispersion and dynamical susceptibility within the continuum theory, it is first useful to 
linearise £ 2 SL [Eq. 33]. This can be accomplished by writing e A (r, r) and e B (r, r) in terms of 8 scalar fields, 

/1 - \ m+ V’a ) 2 + (Vi+r 3 ? + m +^i ) 2 + (v>2 b +^ b ) 2 ] \ 

e A = iR(#h) • V’l + #2 + ^3 + #4 

\ V’l + #2 + V>3 + #4 / 

(\-\ [(^ b - v > 3 b ) 2 + «>! - r 3 ? + (n - V'D 2 + (^ b - ^) 2 ] \ 

e B = in(-0 h ) • ipi + #2 - V’h - ^4 ■ (36) 

\ Vl + i^2 - V’l - #4 / 

The ip* fields are associated with the a bosons in the flavour-wave theory and the ip b fields with the b bosons (see 
Section IID). 

The ip fields are substituted into £ 2 sl [Eq. 33] and terms up to quadratic order are retained. After Fourier transform 
using, 


V>(r,f) 


^3 jduid 2 k e*( k ' r-wt y (k, w), 


(37) 


<S2SL [Eq. A5] can be rewritten as, 


clin 

^2SL — 


l l 


2 (2tt) : 


J dcud 2 k C{ k, u>), 


(38) 


where, 

£( k,w)= * < ( k ,w)-L i (k,w)-tf i (-k,-w). (39) 

i=a,b 

Here, 


^*( k ,w) = (ip\{\t,u),ip l 2 (k,u),ipl(k,bj),ip\(k,u )), 

( Ak ,h + Rk,/i —0 0 > 

Ak ,h Rk,/i 0 0 

0 0 Ak+kM,ft T Rk+kM,/i 

V 0 0 iu> Ak+k M ,/» — Rk+k M ,fo / 


(40) 


L a (k,w) =4 


(41) 
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and, 



(c^ h 

—iijj 

D^h 

0 \ 

4 

iuj 

Ck ,h 

0 

—D^h 


D^h 

0 

C*k+k M ,h 

—iuj 


\ 0 

-D^h 

IUJ 

Ck+k M J / 


(42) 


with, 


A^ h = 4Jn(l + cos 2 2 0 h ) - Jn cos 2 2(9^ k 2 + 2J i2 fc 2 + 4J 22 sin 2 20 h + 2J 22 cos 2 20 h fc 2 


k 2 


B Kh = -4 J n ^1 - — J sin 2 2 0 h + 4 J 22 ^1 - — J sin 2 2 0 h 
^4k+k M ,h = 4Jn sin 2 20 h + Jn cos 2 20 h k 2 + 2 J i2 k 2 + 4J 22 sin 2 2 0 h + 2 J 22 cos 2 26^ k 2 


k 2 


^k+k M J = 4Jn ^1 - 
4Jn ^1 - 

^k+kM,/i = 4Jn ^1 — 

Dkh = —4J 22 (l — 


e 

4 
k 2 

T 

These are just the small k expansions of Eq. 26 and Eq. 28. 

Diagonalisation of Sjsl [Eq. 38] is accomplished using the unitary matrices, 


sin 2 2 0 h + 4 J 22 ^1 - sin2 
cos 2 Oh + 2 Ji 2 /c 2 + 4J 22 
cos 2 Oh + 2Ji 2 /c 2 + 4J 22 
sin 20/i 


U a = 


and, 


ut =vi 


^ ih a (k, cj) itt a (k, cj) 0 0 

-ft a (k, cj) h a (k, cj) 0 0 

0 0 ih a (k + kM,cj) iJ a (k + kM,cj) 

VO 0 -J a (k + k M ,cj) h a (k + kM,cj) 


( iv b (k, uj) iu b (k,uj) iu b (k + kM, uj) iv b (k + kM, uj) \ 

-h b (k, uj) — ft b (k, uj) u b (k + kM,cj) h b (k + kM, co>) 

iJ b (k, uj) — ih b (k, uj) iv b (k + kM, uj) — m b (k + kM, cj) 
y ft b (k, uj) — h b (k, cj) —h b (k + kM,cj) ft b (k + kM, uj) y 


with, 


ft a (k, cj) 


1 






u a (k,w) = 


Bkj, 




+ -BC 


k ,h 


and, 


ft b (k, uj) = 


V2 


\ 




1 + 


C*k ,h C^k+k^jh. 


») +^k, ft 


?i b (k,w) = 


V2 


\ 


C'k,^ — Ck+k M> h 


1 - 


C\c,h C\fL-\-]K.^,h 


0 ) +^M 


In consequence one finds, 


with, 


£(k,w)= ^2 ^*(k,w) • [G*(k,w)] 1 • ^(-k, -u), 

i=a,b 


(43) 


(44) 


(45) 


(46) 


(47) 


(48) 


>T(-k, -w) = uj • <T(-k, -uj) 


(49) 
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Here the diagonal Green’s function matrices are given by, 


[G a (k, w)]' _1 = 4Uj • L a (k, w) • U a = 4 


/ [G a -(k,w)] _1 0 0 0 

0 [G a+ (k, w)] _1 0 0 

0 0 [G a (k + kM,w)] 0 

Vo 0 0 [G a+ (k + kM, w)] ^ 


(50) 


with, 


and, 


[G 3± (k,co)]- 1 =A Kh ±JBl h + 


[G b (k, w)] _1 = 4U* • L b (k,w) • U b = 4 


with, 


(51) 


((k + kM,w)] 1 0 0 0 \ 

0 [G b+ (k + k M , w)] _1 0 0 

0 0 [G b_ (k, w)] -1 0 

\ 0 0 0 [G b+ (k, w)] _1 J 

(52) 


[G b± (k, w)] 1 = -(Ck,h + Gk+k M ,ft) ± -zJ{Ck,h ~ Gk+k M ,/») 2 — 4w(Gk,h - Gk+k M ,ft) + 4 D\ h + 4oj 2 . (53) 


It is useful to calculate the imaginary part of the 
Green’s functions given in Eq. 51 and Eq. 53. One finds, 

^mG a+ (k,w)=0 

QmG a ~ (k, u) = S(ui - u^ h ), (54) 

where, 


,h 




,h ^k ,h 


(55) 


only positive frequency contributions have been retained, 
and, 


^k,/i — y/Alh ^k ,h- 

Similarly, 

^mG b+ (k, uS) =0 

3mG b “(k, w) = S(u ~ u^ h ), 

where, 

T Gk+kM ,h 

y~k,h / _ i 

y (Gk,/i + Gk+k M ,/i ) 2 — h 

and, 

^k,/i = “ Gk+k M ,/i) 

+ 2 \/(Gk,/i + Gk+k M ,/i) 2 — 4-^k,/r 


(56) 


(57) 


(58) 


(59) 


It is clear that ££ h [Eq. 56] and h [Eq. 59] are equiv¬ 
alent at long wavelength to the flavour-wave dispersion 
relations [Eq. 25] and uj^ h [Eq. 27]. 


C. Dynamical spin susceptibility 

The imaginary part of the dynamical spin susceptibil¬ 
ity can be calculated within the linearised theory. Using 
Eq. 12 and Eq. 36 the spin fields can be written as, 

s\(r,t) = — V2(cos 9 h - sm8 h )(ipi + ^ 3 ) 

S y A (r,t) = \/2(cos 9 h + sin6» ?J )('0 b + ^ 4 ) 

S\ (r ,t) = cos 29 h + 2 sin 29 h (V’i + V’i) 

Ss(r,t) = -V2(cos6 h +sm9 h )(ipi - V> b ) 

Sg(r,t) = v/2(cos 9 h - sm9 h )(ip2 ~ ^ 4 ) 

51 (r, t) = cos 2 Oh — 2 sin26 > / l ('0i — / ipl). (60) 

Using the results of Section III B, as well as the definition 
of the spin susceptibility given in Eq. 3, one finds, 

Smx _L (q,w) = 3mx xx (q,w) + 9mx yy (q,w) = 

7 r [(M b (q +q M ,w) cos 9 h -n b (q + q M ,w) sin^) 2 ] 
x 5 b 5(u) - w b ) 

+ 7 T [(n b (q,w) cos 9 h - it b (q, ui) sin^) 2 ] 
x 5 b <5(w-w b +q J 

3mx zz { q,w) = 7 rsin 2 29h [« a (q + qM,w) + n a (q + q M , w )] 2 

X9q +qM ^-w a +qM ). (61) 
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FIG. 7: (Color online). Continuum theory predictions for the imaginary part of the dynamic spin susceptibility of a spin-1, partially- 
polarised, 2-sublattice, spin-nematic state in applied magnetic field. The action c^'sl [Eq. 38] is parametrised using Jn = 1, J 12 = 0.1 and 
J 22 = 2. The external magnetic field is gradually reduced from a) h = 1.5h sa t to f) h = 0. (i) The transverse susceptibility uj) 

[Eq. 61]. Dashed red lines show h at all h and ^q +qM ^ for h < h sat [see Eq. 59]. (ii) The longitudinal susceptibility ^mx zz (q,w) 
[Eq. 61]. Dashed red lines show ^q+q M h [ see Eq. 56]. All predictions have been convoluted with a gaussian to mimic experimental 
resolution. The circuit T-X-M-r in the bond-centred Brillouin zone is shown in Fig. 5. The same linear, normalised colour intensity scale 
is used as in Fig. 6. 
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where either q ^ 0 or q ^ qM- The summed dynamical 
spin susceptibility is plotted in Fig. 7. It can be seen that, 
at long wavelength, this is exactly equivalent to Fig. 6, 
which depicts the flavour-wave prediction for the imag¬ 
inary part of the dynamical spin susceptibility [Eq. 29]. 
Alternatively, this equivalence is clear from directly com¬ 
paring Eq. 29 and Eq. 61. 


D. Hydrodynamic parametrisation 


We now parametrise the linearised field theory c^sl 
[E q. 38] in terms of hydrodynamic parameters. This frees 
the theory from any particular microscopic model, and 
thus allows it to be applied to any partially-polarised, 
2-sublattice AFQ state. We concentrate in particular on 
the case h ~ h sat . 

The action c^sl [Eq. 38] contains all symmetry allowed 
terms for a 2-sublattice, partially-polarised AFQ state on 
the square lattice at a linear level. This is not the case for 
the non-linear action S 2 SL [Eq. A5], which describes the 
long-wavelength fluctuations of the spin-nematic state 
found in [S] [Eq. 4] . For a general 2-sublattice spin- 
nematic state, it may be necessary to include other rel¬ 
evant terms in the non-linear action. However, the only 
effect these will have on the linear theory is to change 
the hydrodynamic parameters. 

In consequence it is possible to write a general, 
hydrodynamically-parametrised, linearised Lagrangian 
for a 2-sublattice, partially-polarised AFQ state, 


^hyd (E, id) q 


Q,Z/ W,Z\ 

Xh 


Q,z 


+ 


S,z S,z / S,z \ 

Xh u k,h\ u - w k \h) 


-|'I' Q ’ z (k, w)| 5 
|§ S ’ z (k, W )| 2 


+ 


uj — uj- 


xy,7r 
k ,h 


|<l ,xy ’ 7r (k, w)| 2 


+ 


UJ — UJ'- 


xy 

k ,h 


■q |# xy,0 (k, w)| 2 


(62) 


This describes 4 modes, and these have all been mapped 
onto k « 0. The dispersion relations of the 4 modes are 
given by, 


"k.k = -J (»? -z ) 2 k2 + fL k4 

-ft=V( A t) 2 +<»t ',) 2k2 

= (WTHOI 2 

= vW°) 2 + W„”) 2k2 - (63) 


on field we write, 

(■ v h ’ Z ) 2 = (! - h/h sat )(v Q ’ z ) 2 
A x y> 7r — a x y>^ ^ 

A xy -° = h, (64) 

while the others are only expected to change weakly with 
varying magnetic field close to h = h S3t . 

The hydrodynamic parameters can be taken from any 
microscopic model supporting a partially polarised, 2- 
sublattice AFQ state. Table I shows the parametrisation 
from ^=i[ S ] [ Eq . 4], 

The [Eq. 63] mode is gapless and associated with 
breaking U(l) symmetry. This is the Goldstone mode 
and describes rotations of the quadrupolar order param¬ 
eter in the plane perpendicular to the applied field. For 
h < h sat the mode has a linear dispersion in the vicinity 
of k = 0. As h h sat from below, (u^’ z ) 2 —^ 0 and there¬ 
fore at h = h sat the mode has a quadratic dispersion, as 
expected for a saturated paramagnet. 

The uj^ z h [Eq. 63] mode is associated with spin fluc¬ 
tuations parallel to the magnetic field. These can be 
thought of as a dynamic spin-density wave. The mode 
is gapped and both the gap, A^ ,z « A^ z , and the ve¬ 
locity, v^ z ~ v^’ z , are only weakly field dependent close 
to h = h sat 5 and thus we approximate them with their 
values at h = h sat . 

The (jJ^h [Eq- 63] mode is associated with spin fluctu¬ 
ations transverse to the field direction and antiparallel on 
the two sublattices. The gap, A^ y,7r , is linearly dependent 
on the field [see Eq. 64], but the velocity u^ y,7r ~ is 
approximately field independent. 

The [Eq. 63] mode is associated with spin fluc¬ 

tuations transverse to the field direction and parallel on 
the two sublattices. The gap, A^ y ’°, is linearly dependent 
on the field [see Eq. 64], but the velocity is 

approximately field independent. 

In order to calculate the imaginary part of the dy¬ 
namic spin susceptibility, it is first necessary to deter¬ 
mine expressions for the spin moments. To do this we 
use ^bbqHS] [Eq. 4] as a guide, and re-express the spin 
fields appearing in Eq. 60 in terms of the fields and hy¬ 
drodynamic parameters appearing in £hyd(k, uj) [Eq. 62], 
making use of the results presented in Section. IIIB. It 
follows that, 


SS x (q, uj) = — sin Oh T xy,7r (—q, —cj) -fcos#/* T xy,0 (—q, — uj) 
^ X (qM + q,cj) = cos 0 h ^r xy ’ 7r (—q, - U ) 

-sin 0 h T xy ’°(-q, -u>) 

ss z {q,tv) = i 2 ^’ z ^ Q,z (-q-~ w ) 

^ z (q M + q,u>) = - sin20 h 4 ,S ’ z (-q,-w), (65) 


For the hydrodynamic parameters that depend strongly 


where q « 0. The imaginary part of the dynamical spin 
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Hydro- 

dynamic 

Spin -1 model 

«b*[s] [Eq- 4 1 

(4 ’ z ) 2 

64Jh(J 22 sin 2 2 6^ + Jn cos 2 2 Oh) 

( A xy, T 

16(Jn — J 22) 2 cos 2 2 9 h 

(A xy ,°)2 

16( Jn + J 22) 2 cos 2 2 6 h 

(vP) 2 

8 (Jn + J 22 X J 11 + 2 J 12 + 2 J 22 ) sin 2 2 6 h 

(*T ) 2 

16 Ji 2 (Jll + J 22 ) + 8 (Jll — J 22 )(Jll — 2J 22 ) 

+ 8 cos 2 2e h [2J 12 (Ju - J 22 ) - (3Ju - J 22 XJ 11 . - 2 J 22 )] 

K y,7r ) 2 

16(J 22 — Jn){Ji2 + J 22 )sin 2 29 h 
+8 Jn (J 22 — Jn) cos 2 20 fl 

K y ’ 0 ) 2 

16(Jn + J 22 ) V 12 + J 22 ) sin 2 2 6 h 
—8 J 11 (J 11 + J 22 ) cos 2 2 6 h 


(Jn + 2 (Ji 2 + J 22 ))/ 3 

x(—Jn + 6 Ji 2 — J 22 + (4Jn + 7J 22 ) cos40/,) 

Q,z 

Xh 

[4(Jn + J 22)]- 1 

S,z 

*h 

KJn ]- 1 

Q,z 

Ph 

2 (Jn + 2 J 12 + 2 J 22 ) sin 2 20h 


TABLE I: Relationship between the hydrodynamic parameters 
appearing in the continuum field theory, Thyd (k, u) [Eq. 62] , and 
the parameters of the microscopic model [S] [Eq. 4] . 


sublattice state realised in Wj Jj 2 [Eq. 1] for J\ < 0, 
J 2 /IJ 1 I > 0.4 and h ~ h sat [see Fig. 1]. In such a state 
the order parameter is bond-centred, but the physical 
spins, the fluctuations of which are measured in inelastic 
neutron scattering experiments, live on the sites of the 
square lattice. 

At low energy, one can view the bond-nematic state in 
terms of an effective spin-1 degree of freedom that lives on 
the bonds 5 . Fluctuations of these spin-1 degrees of free¬ 
dom are described by £h y d(k, uj) [Eq. 62], and, in terms 
of the original spin-1/2 degrees of freedom, correspond to 
changing the mix of triplet states on a bond. However, 
inelastic neutron scattering measures the fluctuations of 
the site-centred spin-1/2, and it is therefore necessary 
to determine the mapping that needs to be applied to 
^mx{ q, a;) [Eq. 66] in order to make experimentally rel¬ 
evant predictions. 

We consider two lattices, the original square lattice 
of the spin-1/2 degrees of freedom, with lattice con¬ 
stant a, and a bond-centred lattice with lattice constant 
b = a/V2 [see Fig. 3]. While in the rest of the article 
the lattice constants have been absorbed into the defi¬ 
nition of the wavevector, rendering it dimensionless, for 
clarity in this section the lattice constants are explicitly 
included in the calculations. In the site-centred lattice 
we label the real-space coordinates by the vector x, and 
in the bond-centred lattice by r. In reciprocal space we 
use p and q. The relations between these coordinates 
are, 


susceptibility is thus given by, 

^rnxbcin, w) = 9mx£(q,w) + 3m^ c (q,w) = 

7T COS 2 Oh S(iO — Uq’h) + Sm 2 0 h S(UJ - LO^h) 

SmXbc(qM +q,w) = 

^ TO x£(qM + q, w) + 3mXbc(qM + q, w) = 

7T COS 2 0 h S(oj - Vq/h) + sin2 “ W q [h) 


3mXbc(q, w) = | ^ ^ S(w - Wq ji) 


7 

q ,h 


SmXb Z c(qM + q, w) = 2ir sin 2 20 h gz gz 


5 ( w - 4h)> 


( 66 ) 


where q ~ 0. We have anticipated the mapping made in 
Section IV and introduced the subscript be which stands 
for bond-centred. 


IV. MAPPING FROM SPIN-1 SITE-BASED 
NEMATIC TO SPIN-1/2 BOND NEMATIC 

We now show how the imaginary part of the long- 
wavelength, dynamical spin susceptibility can be cal¬ 
culated in a bond-nematic phase, focusing on the 2- 


x x — (r x + fy), Xy — 




( r x r y) 5 (67) 


and, 


Px — ^ 2 (^ x + 2y)> Py ~ i/2^ x 


which can be seen from Fig. 5. 

We assume that the spin at a site, S sc , (sc stands for 
site centred) can be written as the sum of the quasi¬ 
spin, S bc , degrees of freedom on the 4 neighbouring bond 
centres, which leads to, 


QSC _ ebe 


?X C • 


where e x = (a, 0) and e y 
transform results in, 


ebe 

' Xi + ir 


(0, a). 


+ S'x'Lfy > (69) 

2 

Taking the Fourier 




Px ~ 


Qx 


V2 


q y 2rmr\ 




X 6 ( Py + 


V2 


q y 2nir\ q x b q y b 

— ± - I ms-ms —— 


a / 


COS — COS -y,S q 


be 


(70) 


where m and n are integers. This allows the imaginary 
part of the dynamical spin susceptibility to be calculated 
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FIG. 8: (Color online). Predictions for the dynamic spin susceptibility of a spin-1, partially-polarised, 2-sublattice, spin-nematic state in 
applied magnetic field (Fig. 7) mapped onto the site-centred lattice The relationship between the bond- and site-centred lattices is shown 
in Fig. 5, and the mapping of the dynamical susceptibility is performed using Eq. 71. The external magnetic field is gradually reduced 
from a) h = 1.5h sa t to f) h = 0. (i) The transverse susceptibility Srnxi(q, w) [Eq. 71]. Dashed red lines show h and ^q +qM ^ [ see 
Eq. 59]. (ii) The longitudinal susceptibility S7nx zz (q, cj) [Eq. 71]. Dashed red lines show h and ^q+ qM h [ see 56]. All predictions 

have been convoluted with a gaussian to mimic experimental resolution. The circuit T-X-M-r in the site-centred Brillouin zone is shown 
in Fig. 5. The same linear, normalised colour intensity scale is used as in Fig. 6. 
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in the site-centred coordinate system as, 


9m Xs T(p,^) = 
2 


COS 


■ sin 


(Px-Py)a' 

^ 2 

’fax +Py)a 

4 

cos 

4 

*5" 

X 

1 

< 

77 

•2 

’(Px +Py)a 

4 

sin 

4 


Smx b“(P»w) 

^Xbc“(PM +P,w) 

(71)' 


where (p x + p y ) 2 /2 + (p x — p y ) 2 /2 = p 2 has been used, 
p « 0 and 9mXb c a is given in Eq. 66. 

The result of mapping the spin-1 dynamical suscepti¬ 
bility predictions shown in Fig. 7 onto the site-centred 
lattice using Eq. 71 is shown in Fig. 8. It can be 
seen that the mapping from the bond-centred to site- 
centred lattice results in all the low-energy modes ap¬ 
pearing at the T point, which is expected, since the 
2-sublattice AFQ state does not break the transla¬ 
tional symmetry of the site-centred lattice. For q « 0, 
one has cos 2 [(g x — q y )a/4\ ~ 1, cos 2 [(g x + g y )a/4] ~ 1, 
sin 2 [(g x — q y )a/4] « 0 and sin 2 [(g x + q y )a/ 4] « 0 and 
therefore the contribution from (q, oj ) dominates 

over the contribution from Qm + q, oS). In con¬ 

sequence, when making experimental predictions in Sec¬ 
tion VI below, we concentrate on the modes l and 

[Eq. 63], since these will determine the dominant 
experimental signatures. 


V. MICROSCOPIC THEORY OF A SPIN-1/2 
BOND NEMATIC 


quadrupolar ordering of the spin degrees of freedom just 
below the saturation magnetic field, and hence the cre¬ 
ation of a spin-nematic state. We show below that this is 
the case for / Hj i "lY 2 2 [Eq. 1] over a wide range of parame¬ 
ters. 

In the saturated paramagnet it is useful to use the 
hardcore boson representation, in which spin operators 
are replaced according to, 

Sf =-1/2 +a\ai , s+=a\, Sf s= a t . (72) 

Thus [Eh- 1] can be rewritten as, 

= -sl = “ MM] a q«q 

q 

+ 2JV (73) 

V q,k,k' 

where, 

cj(q) = e(q) - e min , [i(h) = h cl - h, 

h c i = e(0) - 6 min , V q = 2(e(q) + U), (74) 

U —y oo is the on-site repulsion and, 

e(q) = J i (cos q x + cosg y ) + 2 J 2 cosq x cosq y . (75) 

In the parameter range — 2 < J 1 /J 2 < 2 with J 2 > 0, 
6 m in = e(Q±) = —2J 2 , where Q + = (7r,0) and 
Q_ = (0,7r). The saturation field for 1-magnon 
condensation is thus given by, 


In order to make quantitative experimental predictions 
for materials it is necessary to determine the hydrody¬ 
namic parameters appearing in £h y d(k, ui) [Eq. 62]. Here 

we consider [Eq. 1] and microscopically calculate 

the hydrodynamic parameters in the vicinity of h = h sat . 
In Section VI these will be fed into £h y d(k, uS) [Eq. 62], 
allowing quantitative predictions to be made for inelastic 
neutron scattering experiments. 

All magnets have a saturation magnetic field, above 
which the spins are aligned parallel to the field direction. 
For systems with an antiferromagnetic ground state at 
h = 0, the simplest way to connect the zero-field and 
high-field states is via a canting of the ordered moment 
towards the field direction. As field is reduced from above 
saturation, there is a phase transition from the fully- 
polarised state to the canted antiferromagnet at the sat¬ 
uration field. This can be understood in terms of the 
condensation of magnons out of the fully-saturated “vac¬ 
uum” state 53,54 . 

In frustrated magnets the condensation of single 
magnons may not be the first instability of the fully- 
saturated state as magnetic field is lowered. One possi¬ 
bility is that the 1-magnon instability is preceded by a 
2-magnon instability, in which bound pairs of magnons 
condense 5 . If this occurs, then this corresponds to a 


h cl =2Ji+4J 2 . (76) 


For 11 (h) < 0 (ft, > h c 1 ) the system is fully-polarised 
unless the 1-magnon instability is preceded by another 
instability. We denote the fully-saturated state by |fi), 
where ai\Q) = 0. Since the interaction term in Eq. 73 
is normal ordered, there is no self-energy term in the 1- 
magnon dispersion relation, and it is exactly given by, 


= w ( k ) - mW- 


(77) 


Expanding o;(k) near the minima Q± leads to, 


j(k x , ky + 7r) — -Tjy H-yyr + 0(k 4 ) 

2 m[ L) 2^2 


(78) 


where, 


m 


(i) 


A 1 ) 


1 


1 -Ji + 2J 2 ’ 2 


m 0 = 


J\ + 2 J 2 


(79) 


and, 


u{k x + 7 r, k y ) = H-jry + 0(k 4 ). 


2 rri 2 2 m\ 


(80) 
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r(A, K;p, p ) = 


K/2+p K/2+p' K/2+p 

+ 


K/2+p" K/2+p' 


r(A, K;p, p ) 


K/2-p K/2-p' K/2-p 


K/2-p" K/2-p' 


FIG. 9: (Color online). Diagrammatic calculation of 2-magnon 
binding energy for [Eq. 1]. The exact scattering amplitude 

T(A,K;p,p / ) [Eq. 82] is calculated from a ladder diagram. K is 
the centre of mass momenta and A the binding energy of magnon 
pairs. 


In the fully-saturated state it is very simple to calculate 
the imaginary part of the dynamical spin susceptibility. 
It is exactly given by, 

( w > k ) = k ) + 3raXsat(w, k ) = 

3m Xsat( w > k ) = °- ( 81 ) 

This describes a sharp band of 1-magnon excitations, 
with equal weight at all wavevectors. 

We next consider the condensation of bound pairs of 
magnons out of the fully-saturated state. The 2-particle 
Green’s function can be calculated exactly using the lad¬ 
der diagram shown in Fig. 9. The scattering amplitude 
is given by 19,24,53-55 , 

r(A, K; p, p') = Vp’-p + V-p'-p 

_1 [ dV r(A, K; p, p")[yp/_p// + V_p'_ p »] 

2 J (2t r) 2 ui( K/2 + p") + w(K/2 - p") + A - i0+ ’ 

(82) 



JM\ 

FIG. 10: (Color online). First instability of the fully-polarised 
state on lowering magnetic field. The critical field for one magnon 
condensation [h c i, Eq. 76] is shown as a red, dashed line and the 
critical field for condensation of bound-magnon pairs [h c 2 , Eq. 83] 
as a blue line. For 0.4 < | J 1 /J 2 I the 2-magnon instability precedes 
the 1-magnon instability (h c 2 > h c 1 ). 



where A and K are respectively the binding energy and 
the center-of-mass momentum of the bound state. This 
integral equation is exactly soluble 53,54 . 

Divergence of r(A,K;p,p / ) [Eq. 82] implies the ex¬ 
istence of a stable bound state below the 2-magnon 
continuum 19,56 . The K-dependent binding energy is de¬ 
noted as A#(K). The wavefunction of the bound state 
can be determined from the residue of the divergence of 
r(A, K; p, p') [Eq. 82], and this is considered in more de¬ 
tail below. When lowering magnetic field, if the bound- 
magnon gap closes before the 1-magnon gap then the 
bound state condenses and the spin-nematic state ap¬ 
pears. The critical field for 2-magnon condensation is 
given by, 

h C 2 = h c i + A m /2, (83) 

where A m is the maximum value of the binding energy. 
For 0.4 < IJ 1 /J 2 I < 5 there is a stable bound state at 
K = (0,0), and this is the leading instability on low¬ 
ering magnetic field [see Fig. 10]. The binding energy, 
A m = A#(K = (0,0)), as a function of J 2 /IJ 1 I is shown 
in Fig. 11. 


FIG. 11: (Color online). The 2-magnon binding energy in the 
S=1 /2 

fully-saturated state of [Eq. 1]. The binding energy, 

Am = A j g(K = (0,0)), is calculated from r(A,K;p,p') [Eq. 82, 
Fig. 9]. It determines the saturation field for 2-magnon condensa¬ 
tion, h c 2 = h c 1 + A m /2. For A m > 0 [see Fig. 10], h c 2 > h c 1 and 
the 2-magnon condensate occurs at higher field than the 1-magnon 
condensate. 


In the fully-polarised phase with h > h c 2 , the bound- 
magnon dispersion relation is given by, 

<4T = A m -A B (K)- M2 (/ i ) * Afi- M2 (ft) + 0(K 4 ), 

(84) 

where, 

112 (h) = 2(h c2 - h). (85) 

This dispersion, [Eq. 84], is shown in Fig. 12 in 

relation to the 2-magnon continuum. 

Slightly below the saturation field, for small 112 (h) > 0, 
one can view the system as a dilute gas of bound 
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FIG. 12: (Color online). 2-magnon dispersion of 4^ 

[Eq. 1] at saturation. The 2-magnon continuum, defined by 
lj( K/2 + p) + ou(K/2 — p) + A m , is shown as a blue region. The 
dispersion of the 2-magnon bound state, [Eq. 84], is shown 

as a white dashed line for the region of K in which it lies below the 
2-magnon continuum, (a) shows the full extent of the continuum, 
while (b) shows a detailed view of the low-energy spectrum. The 
exchange parameters are taken as J\ = — 3.6K and J 2 = 3.2K and 
the field as h = h c 2 [Eq. 83] . 


magnons. We consider an effective Hamiltonian, 
'H e s= ^2 { 2 m ( 2 ) - j b K b K 

|K|<A k ' 

p(2) _^ 

+ TV 52 ^Kl+q & K 2 -q &K i &K 2 4 

Ki,K 2> q 

where the bound-state creation operator is, 

= XK (p) a K/2+p a K /2-P 7 


( 86 ) 


(87) 


where 6(x) — jj f (PKb-^e 11 ^ ^. Substituting, 

6(x) = y/p 2 +Sp 2 (x) e* e(x) , (91) 

into <S e ff [Eq. 86] and integrating out the high-energy 
mode Sp 2 , results in, 

5 eff = Jdx 2 ^(d t 6{x)) 2 - ^g_(V0(x)) 2 . (92) 

Hence, the disperison of the Goldstone mode at small K 
within the dilute limit is given by 57 , 


2—m 

K, 


r( 2 y 2 K 2 

2 m( 2 ) 



|K|. 


(93) 


M r <2) /\Ji\=2i 



FIG. 13: (Color online). Magnetisation of the spin-nematic state 
found in 4 ^ Ji _j 2 [Eq. 1] just below the saturation field as a function 
of magnetic field h. M(h) [Eq. 94] is plotted at several values of the 
effective interaction O 2 ). Comparison of dM(h)/dh at h —* h sa t 
[Eq. 95] with experimental data can be used to estimate G 2 ) for a 
given material. 


is the renormalized interaction between the low- 
energy bound magnons 24 , A is a momentum cutoff, and 
the dots represent higher-order interaction terms that can 
be neglected in the dilute limit. The effective free energy 
of the condensed phase is given by, 

^ -q r(2 V! -P2{h)p 2 , (88) 


In 2d, the parameter T^ 2 ) tends to zero approaching 
the saturation field. However, even a very small inter¬ 
layer coupling will cutoff this approach to zero, and T^ 2 ^ 
remains finite at h = h c 2 58 . A good way to estimate 
F( 2 ) for a given material is from the gradient of the ex¬ 
perimental magnetisation curve close to saturation. One 
finds theoretically for the magnetisation, 


where (6 k=o) 
results in, 


s/p 2 e ie °. Minimizing this free energy, 


_ 2 ii 2 (h) 

p2 - r( 2 ) 


(89) 


M(/. ~ M * j) I><* > = i- 16( rm (94) 

and therefore, 


Next we consider the Goldstone mode of the spin- 
nematic state in the dilute limit at zero temperature. 
Taking the cutoff A —)► oo, the effective action is given 
by, 


/■ 


<S e ff = / d 2 x -b(x)*d t b(x) - b(x)d t b(x)*) 
o r( 2 ) 

+M2 (h)lb(x)l 2 -—lb(x)l 4 


|VKx)| 2 

2m( 2 ) 

(90) 


dM 

~dh 


h — /l S at 


16 

fo) - 


(95) 


Magnetisation curves for several values of T^ 2 ^ are shown 
in Fig. 13. 

Finally, we will consider the bound-state wavefunction 
Xk(p) [Eq. 87] and demonstrate that this describes the 
same 2-sublattice bond-nematic state considered in Sec¬ 
tion II and Section III. 
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From taking the residue of the divergent scattering am¬ 
plitude r(A,K;p,p') [Eq. 82], one finds, 


Xk=o(p) 


Xo(p) oc 


COS p x — COS Py 
2cj(p) + Am 


(96) 


where the minimum of cj(p) [Eq. 74] is zero. The wave- 
function is normalised by requiring ^ Ixo(p)| 2 = 1, 
where the summation over p is taken in half of the Bril- 
louin zone as p and — p gives the same bound state. This 
wavefunction has d-wave symmetry, since interchange of 
the coordinates x and y leads to Xo(p) — Xo(p) while 
K = 0 —y 0. 

The bound-magnon condensed phase is described by 
the coherent state 16 , 


|Nem) = Ciexp(^y]xo(p)aJ ) a-p)|fi>, (97) 

P 


where C\ = n pV /l - |</>Xo(p)| 2 and 0 = ^/p 2 e z9 ° . In 
this phase we consider the expectation value of the bond- 
operator, 

Q ( r ) = ^1 ^l+r = 2 (^M+r “ Cp+r) “ ^lfl+r 
= Q{l7r ~ ^M+r’ ( 98 ) 

where is defined in Eq. 2 and 

Q*j~ y =1/2 (Q™ — Q™) is defined in analogy with 
Eq. 6 for the spin-1 theory. 2-sublattice AFQ order of 
the type considered throughout this article corresponds 
to a non-zero expectation value (Q (r)) on nearest- 
neighbour bonds, with a change of sign between vertical 
and horizontal bonds. Calculating the expectation value 
of this bond operator with respect to the condensed 
phase, |Nem) [Eq. 97] gives, 


{Q (l*))nem — ($\ ^i_|_ r )nem — ^ ^ ((2p(2_p) n em exp(ip • I*) 

P 


_ ^Xo(p) 

„ 1 “ l0Xo(p)| 2 


exp(ip • r) 


0^ X °( p ) ex P(*P-r) + O(0 2 ). (99) 


The permutation r = (r x ,r y ) —» (r y ,r x ) shows the d- 
wave symmetry of this bond operator, 


{Q (?xi ^y))nem — {Q nem 5 (100) 


which, on nearest-neighbour bonds, exactly corresponds 
to the 2-sublattice AFQ phase shown in Fig. 1 and Fig. 3. 

Next, we consider the asymptotic behaviour of 
(Q __ (r))nem [Eq. 99] for the low-density case </> <C 1. 
For large r, the oscillation of the wave function is fast, 
and the dominant integrand comes from p close to (0, 0), 
(0,7r) and (7r,0). Considering r = (r, 0) for simplicity, 
one can show for r 1, 


(Q (r, 0))nem 


C2 


exp(— r/go) 

\X7Co 


( 101 ) 


where £ 0 « l/y2m^A m is a measure of the size of the 
bound state and o 2 is a constant. The larger the value of 
the binding energy A m [see Fig. 11], the more localised 
the bound state. 

As an example, we take J 2 /IJ 1 I = 0.9 and show that 
the nearest-neighbour bonds gives the dominant contri¬ 
bution to (Q (r)) [Eq. 99]. One finds, 

(Q—(2,0)) nem /(Q—(l,0)) nem = 0.15 
(Q— (3,0))nem/(Q”(l,0))nem = 0.33 
{Q— (4,0)) nem /(Q— (l,0)) nem = 0.10 
{Q— (0,0)) nem = (Q“(l, l)) ne m = 0 (102) 

to first order in </>. In Section II, we have in ef¬ 
fect considered the nearest-neighbour order parameter 
(Q -- (l, 0))nem> and the above analysis shows that this 
approximation becomes better the larger the value of 
A m . In the limit where A m —»• oo the wavefunction 
Xk=o(p) [Eq. 96] is only non-zero on nearest-neighbour 
bonds to first order in 0, and we write, 

Xn.n.(P) = \fj^( COS Px ~ cos Py)- (103) 

In consequence the only non-zero bond-operator expec¬ 
tation values are, 

(Q—(l,0)) nem = -(Q“(0, l)) ne m = 0- (104) 

While the mapping to the spin-1 model [Section II] 
needs a nearest-neighbour bond order, this is not the case 
for the continuum theory [Section III]. The sole require¬ 
ment is K <C l/£o> and therefore the continuum theory 
is valid even for small A m . 

In summary, we have considered in this section the 
ground state and low-energy excitation spectrum of 
Wj"! 1 / 2 [Eq. 1] for magnetic fields close to saturation. 
At the saturation field, h sat = h c 2 , bound-magnon pairs 
condense to form a spin-nematic state. In the limit 
that the density of bound pairs is dilute, the excitation 
spectrum can be calculated from the microscopic model. 
This allows the hydrodynamic parameters appearing in 
£hyd(k, cj) [Eq. 62] to be determined, and we show an 
explicit example of this in the following section. 


VI. PREDICTIONS FOR INELASTIC 
NEUTRON SCATTERING EXPERIMENTS IN A 
SPIN-1/2 BOND NEMATIC 


Spin-nematic order does not break time-reversal sym¬ 
metry, and therefore does not produce an internal mag¬ 
netic field. This makes the spin-nematic state essentially 
invisible to most common probes of magnetism, such as 
elastic scattering of neutrons, Knight shift of the NMR 
spectra and the asymmetry of oscillations in jusr spectra 3 . 
However, excitations of the quadrupolar order parameter 
mix a spin-dipole component into the wavefunction, and 
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FIG. 14: (Color online). Predictions for inelastic neutron scatter- 
ing from a system described by "Hj~_ j 2 [Eq. 1] above the saturation 
field h sa t = h C 2 [Eq. 83]. The imaginary part of the dynamic spin 
susceptibility, Smx s ^ t (w,q) [Eq. 81], is calculated exactly using 
the parameters J\ = — 3.6K and J 2 = 3.2K, which are believed 
to describe BaCdV 0 (P 04 ) 2 25 . The 1-magnon dispersion, 

[Eq. 77] , is shown by a red dashed line and the dispersion of the 2- 
magnon bound state, [Eq. 84], by a white dashed line, a) At 

h = 1.23hsat the dispersion of the 2-magnon bound state is gapped 
for all q. b) At h = h sa t the gap to the 2-magnon bound state closes 
at q = (0, 0) while 1-magnon excitations remains gapped for all q. 
The circuit r-X-M-T in the site-centred Brillouin zone is shown in 
Fig. 5. The same linear, normalised colour intensity scale is used 
as in Fig. 6. 


this can, in principle, be detected by dynamic probes of 
magnetism. 

Here we make predictions that demonstrate how a 
spin-nematic state could be identified via inelastic scat¬ 
tering of neutrons. We consider in particular the 
class of materials that are well described by 
[Eq. 1]. For definiteness we show predictions relevant 
to BaCdV 0 (P 04 ) 2 , where fits to magnetic susceptibility 
give the exchange parameters J\ = — 3.6K and J 2 = 3.2K 
[25]. We take a square lattice with lattice constant 
a = 4.5A. The saturation field of BaCdV 0 (P 04)2 has 
been measured to be ft sat ~ 4.2T, which is easily achiev¬ 
able in a neutron scattering experiment. 

In order to make quantitative predictions for experi¬ 
ment it is necessary to determine the hydrodynamic pa¬ 
rameters appearing in Smx s a c a (p, uj) [Eq. 71]. All mo¬ 
mentum transfers are now relabelled as q, since this is 
commonly used. We concentrate on the 1-magnon mode 
with dispersion [Eq. 63] and the 2-magnon mode 

with dispersion [Eq. 63], as these will be the most 
experimentally visible. 

Comparing [Eq. 63] to [Eq. 84 and Eq. 93] 

and using the coefficient of Seff I E( h 92 1 to Set Ph ,Z > ° ne 


finds for the 2-magnon mode, 
1 


a h S at — 


( X ?' 1 ) 2 = 


2 ra( 2 ) 

fi 2 (h) 2 ft s , 


rn 


( 2 ) 


rn 


( 2 ) 


1 


h & 


q,z _ ft2(ft) _ 2ft sat 

^ h r ( 2 )m( 2 ) r( 2 )m( 2 ) 


1 - 


Q, z 

xi' 1 = 


1 


( x?’ ! ) 2 r< a) ' 


(105) 


where ft sat = h c2 . For the 1-magnon mode one can com¬ 
pare [Eq. 63] to [Eq. 77] to find, 


A7’° = h 


«:. 0 ) 2 = -- 


Ji + 2J 2 


(106) 


For J 2 /\Ji\ = 3.2/3.6 « 0.89, one can use the results of 
Section V to show, 


rn = 1 ' 56, rn = °' 196, 

Mil Mil 

TyT = 1-65, = 1.27. 

Mil 


(107) 


Furthermore, from fitting the experimentally mea¬ 
sured magnetisation curve for BaCdV 0 (P 04 ) 2 25 with 
M{h « h sat ) [Eq. 94], one can estimate, 


p( 2 ) 

n\ 


10 . 


(108) 


For h > h sat it is possible to exactly calculate both 
9mX s “?(w,q) [see Eq. 81] and the 2-magnon dispersion 
relation, [Eq. 84]. In Fig. 14 we show predic¬ 

tions for inelastic neutron scattering at h = 1.23h sat and 
h = h sat . The only signal is a sharp and uniformly intense 
band of 1-magnon excitations. Also shown in Fig. 14 is 
the 2-magnon dispersion, which is gapped for h > h sat , 
and softens at q = 0 for h = h sat , but is invisible to in¬ 
elastic neutron scattering experiments. In Fig. 15 these 
predictions are integrated over 47r of solid angle in or¬ 
der to mimic inelastic neutron scattering from a powder 
sample. 

For h < h sat we calculate predictions for inelastic neu¬ 
tron scattering using Eq. 71. These predictions are shown 
in Fig. 16 at a range of magnetic field values and in the 
vicinity of q = 0. 

As field is reduced below h = h sat intensity appears 
in the Goldstone mode excitation, [Eq. 63], for 

q 7 ^ 0. For small q and fixed ft, the intensity of 
this mode grows linearly with q, as can be seen from 
Eq. 71. For fixed, small q and small ft sat — ft the intensity 
grows as V^sat — ft- The velocity of this mode goes as 
v h ’ z V^sat ~ ft [Eq. 105], and therefore the dispersion 
becomes steeper as field is reduced. 
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FIG. 15: (Color online). Predictions for angle-integrated inelastic neutron scattering experiments on a spin-1/2 frustrated ferromagnet 
with incipient spin-nematic order, at the saturation field h = h sa t- The imaginary part of the dynamic spin susceptibility [Eq. 81] has been 
integrated over 47r solid angle in order to mimic a powder sample. The Hamiltonian considered is 'Hj/Zj 2 [Eq. 1], and the parameters 
Ji = —3.6K and J 2 = 3.2K are believed to describe BaCdV 0 (P 04 ) 2 25 - The white dashed line shows the dispersion of the 2-magnon 
continuum, [Eq. 84], which becomes gapless at the saturation magnetic field h = h sa t- At this value of magnetic field, spectral 

weight resides in the 1-magnon excitation, which is gapped for all momentum transfers, q. The same linear, normalised colour intensity 
scale is used as in Fig. 6. 


The 1-magnon excitation has a gap A^ y, ° = h that 
slowly reduces with lowering field. To a first approxi¬ 
mation the velocity is constant, and the intensity in this 
mode is uniform over q and does not vary with h. 

The relative intensity of the two modes can be esti¬ 
mated from, 


jQ.zfq, h) _ q 2 

tcy,o(q,/i) ~ 


(109) 


where 7q >z is the intensity of the Goldstone mode and 
I xy ,o the intensity of the gapped 1-magnon mode. 

In Fig. 2 we show predictions for inelastic neutron scat¬ 
tering from a powder sample, calculated by taking the 
predictions shown in Fig. 16 and averaging over 47r of 
solid angle. We show predictions for h = 0.8h sat . Also 
shown is a constant-q cut at q = 0.18A -1 , showing the 
relative intensity of the two modes. At these values of q 
and h, the peak intensity in the Goldstone mode is about 
3% of the peak intensity in the 1-magnon mode. 

There will also be a small contribution to the inelastic 
scattering from the 2-magnon continuum [see Fig. 12]. 
This is spread over a large region of q and uj space, and 
at leading order the intensity grows as h sat — h as the field 
is reduced. Thus the contribution to the scattering will 
be considerably smaller than from the Goldstone mode 
excitation, which is sharp and has an intensity growing 
as \/hsat ~ h, and can be safely ignored. 

While the predictions shown in Fig. 14, Fig. 16 and 
Fig. 2 are specific to BaCdV 0 (P 04 ) 2 , a very similar anal¬ 
ysis can be made for any compound described by 
[Eq. 1]. In fact, £h y d(q> Co ’) [Eq. 62] can be applied to 
any system with a partially polarised, 2-sublattice spin- 
nematic order parameter. 


VII. DISCUSSION AND CONCLUSIONS 


In this article we have explored how inelastic neutron 
scattering can be used to probe for the existence of a 
spin-nematic state in applied magnetic field. Following 
the philosophy detailed in Refs. [21-23], we suggest that a 
good way to recognise this state experimentally is via the 
excitation spectrum, since the ground state is essentially 
invisible to common probes of magnetism. To this end, 
we have developed a general theory of the magnetic ex¬ 
citations of a two-sublattice, bond-centered spin-nematic 
state in applied magnetic field. We parameterise this the¬ 
ory from the microscopic model believed to describe the 
spin-1/2 frustrated-ferromagnet BaCdV 0 (P 04)2 [25], a 
promising candidate for spin-nematic order, and make 
predictions for inelastic scattering of neutrons from this 
material. We also introduce an effective spin-1 model 
supporting the same form of spin-nematic order, and use 
it to explore the evolution of magnetic excitations for 
a wide range of magnetic field. The main experimental 
predictions of this article are summarised in Fig. 2 and 
Fig. 16. 

The starting point was to first derive a phenomeno¬ 
logical theory of a 2-sublattice, partially-polarised spin- 
nematic state. This involved constructing a spin-1 
bilinear-biquadratic model [77[S], Eq. 4] with both 
the desired symmetries and a spin-nematic ground state. 
We note that if spin-1 compounds with large biquadratic 
coupling can be synthesised, this model may become ex¬ 
perimentally relevant in its own right. Alternatively, 
it could be realised in molecular condensates of cold 
atoms. In this article, the spin-1 model served as a 
guide to the derivation of a continuum field theory de¬ 
scription of the long-wavelength excitations of the spin- 
nematic state [£h y d(q, Eq. 62]. We found four low- 






23 


(i) Qmx ± (q, w) 


(ii) Smx zz (q,w) 



04 ) 
0.00 


0.05 


0.10 0.15 

q ( 1 /A) 


0.20 


0,0 

0.00 


0.05 


0.10 0.15 

q ( 1 /A) 


0.20 


0,0 

0.00 


0.05 


0.10 0.15 

q ( 1 /A) 


0.20 



0,0 

0.00 


0.05 


0.10 0.15 

q ( 1 /A) 


0.20 


0 , 0 ' 

0.00 


0.05 


0.10 0.15 

q ( 1 /A) 


0.20 



b) / 7 = 0 . 9 / 7 sat 


0 , 0 ' 

0.00 


0.05 


0.10 0.15 

q ( 1 /A) 


0.20 



c) h= 0 . 8 h sat 


o,p' 

o.oo 


0.05 


0.10 0.15 

q ( 1 /A) 


0.20 


0.61 

0.5 

0.4 

0.3 

0 . 2 | 

0.1 


o,o' 

0.00 


d) h= 0 . 7 h sar 


0.05 


0.10 0.15 

q (1/A) 


0.20 


1. 


2.5x1 O’ 1 


6.3 xIO" 2 


1.6 x 10 -2 


4. xIO -3 


l.xlO' 3 


FIG. 16: (Color online). Predictions for inelastic neutron scattering from a single-crystal sample of a spin-1/2 frustrated ferromagnet 
exhibiting two-sublattice, bond-centered spin-nematic order in applied magnetic field, (a)—(d) Predictions at small q, for magnetic field 
h ranging from the saturation value h sa t to h — 0.7h sa t- (i) The majority of spectral weight resides in a gapped spin-wave mode, visible in 
the transverse part of the dynamical susceptibility <*>)• (ii) The onset of spin-nematic order for h < h sat is heralded by the emergence 

of a ghostly, linearly-dispersing Goldstone mode in the longitudinal susceptibility Predictions for X° i ^( c L^ UJ ) were calculated for 

a material described by the spin-1/2 J 1 —J 2 model [Eq. 1], as described in Section VI of this article, for parameters J\ = —3.6 K 

and J 2 = 3.2 K relevant to BaCdV 0 (P 04 ) 2 25 . All predictions have been convoluted with a gaussian of standard deviation 0.006 meV to 
mimic experimental resolution. Equivalent results for a powder sample are shown in Fig. 2. 
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energy modes, including a Goldstone mode associated 
with the breaking of U(l) symmetry and predominantly 
describing quadrupole fluctuations of the order parame¬ 
ter, and three gapped modes describing mixed spin and 
quadrupole fluctuations. 

One of the most experimentally promising places to 
search for the spin-nematic state is in square lattice, spin- 
1/2 frustrated ferromagnets described by Wjl 1 / 2 [Eq. 1]. 
Close to saturation the spin-nematic state is expected to 
be realised for a large range of J 2 /J 1 values. In spin-1/2 
nematic states the quadrupolar order parameter lives on 
the bonds of the lattice and not on the sites. The con¬ 
tinuum theory thus describes effective bond-centred spin 
fluctuations, and in order to accurately describe exper¬ 
iments, it was necessary to make a mapping onto the 
site-centred lattice. After this procedure it became clear 
that the experimental response is dominated by only 2 
modes, the Goldstone mode and a gapped mode that can 
be identified with single-magnon condensation out of the 
fully-saturated state. 

The continuum theory contains a number of hydrody¬ 
namic parameters that have to be calculated from mi¬ 
croscopic considerations. In order to do this we consid¬ 
ered Hjl 1 / 2 [Eq. 1 ] close to the saturation magnetic field, 
and used exact diagrammatic calculations to determine 
the saturation magnetic field, assuming condensation 
of first single magnons and then bound-magnon pairs. 
For a wide range of J 2 /J 1 values the bound-magnon 
pairs condense at higher field than single-magnon exci¬ 
tations, forming a spin-nematic state. Using similar dia¬ 
grammatic calculations, the continuum theory was fully 
parametrised. 

This allowed us to make predictions for inelastic neu¬ 
tron scattering experiments, focusing on the material 
BaCdV 0 (P 04 ) 2 - In Fig. 2 and Fig. 16 we show predic¬ 
tions for small momentum transfers q and fields h < h sat , 
where the dominant feature in the spectrum is the gapped 
1-magnon excitation band. However, we showed that 
there is also spectral weight in the Goldstone mode ex¬ 
citation and this grows with increasing q and decreasing 
h. Experimental detection of this excitation would be 
strong evidence for the existence of a spin-nematic state. 

The predictions we make for BaCdV 0 (P 04)2 are 
quantitative, and this allows one to determine whether 
it is realistic to expect to see the Goldstone mode ex¬ 
perimentally. The saturation field has been measured as 
h sat ~ 4.2T, and therefore the full field range 0 < h < h sat 
can be accessed in inelastic neutron scattering exper¬ 
iments. Assuming this material is well described by 
[Eq. 1 ], one would expect a phase transition at 
intermediate field between the spin-nematic state (higher 
field) and a canted antiferromagnet (lower field). Since 
the frustration parameter J 2 /I Ji| ~ 0.9 is relatively close 
to the highly frustrated point J 2 /IJ 1 I ~ 0.5, one would 
expect that the spin-nematic state would be the ground 
state over a sizeable field range. The lower the mag¬ 
netic field in which the spin-nematic can be measured, 


the more intense the scattering from the Goldstone mode. 
Taking a relatively conservative value of h = 0.8h sat [see 
Fig. 2] one can see that an energy resolution of better 
than about O.lmeV is needed at a momentum transfer 
of about 0.15 — 0.2 A -1 - achievable values in neutron 
scattering experiments. The intensity of the Goldstone 
mode at these values of field and momentum transfer 
is expected to be about 3% of the 1-magnon excitation. 
This is comparable to measuring an ordered moment of 
about 0.2/iB - and this is possible experimentally. 

Finally we would like to emphasise the generality of 
these results for understanding spin-nematic phases, and 
end with the hope that the ghostly Goldstone mode will 
be revealed through the machine of inelastic neutron 
scattering. 
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Appendix A: Continuum theory of the 2-sublattice 
antiferroquadrupolar state at h=0 

At h = 0 the continuum field theory developed in Sec¬ 
tion III for the 2 -sublattice antiferroquadrupolar (AFQ) 
spin-nematic state can be considerably simplified. This 
allows comparison to lattice gauge theory calculations 
used to understand the h = 0 spin-nematic state in 

UA 2 [Eq- I] 33 " 35 - 

One simplification arises from the increased symmetry 
at h = 0 . In the case h 7 ^ 0 the symmetry of <S 2 Si_ [Eq. A5] 
is U(l), and breaking this symmetry via the formation 
of an AFQ state results in a single Goldstone mode, as 
shown in Fig. 7b-e. For h — 0, the symmetry of the action 
is increased to SU(2). As a result, there are three Gold¬ 
stone mode excitations, as can be seen in Fig. 7f) There 
is a degenerate pair of Goldstone modes associated with 
rotations of the quadrupolar order parameter out of the 
ordering plane (Fig. 7(i)f), and a third associated with 
rotations within the ordering plane (Fig. 7(ii)f). There 
is also a gapped mode that can be interpreted as a dy¬ 
namical spin-density wave excitation. 

A second simplification arises from a natural division 
at h = 0 of conjugate pairs of fields into those asso¬ 
ciated with high- and low-energy fluctuations. This is 
not the case as h h sat . At h == 0 , fields associ¬ 
ated with high-energy fluctuations can be eliminated by 
a Gaussian integral. The resulting continuum theory is 
based on an SU(3) generalisation of the non-linear sigma 
model (nlcrm). A similar treatment of the 3-sublattice 
AFQ state on the triangular lattice was presented in 
Refs. [21,23], and the development in this Appendix 
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FIG. 17: (Color online). 2-sublattice antiferroquadrupolar (AFQ) 
spin-nematic state on the site- and bond-centred lattices at h = 0. 
Green spheres represent spin-1/2 degrees of freedom at the vertices 
of a square lattice (site-centred lattice). The bond-centred, nematic 
order parameter is represented by red cylinders, and the associated 
distribution of spin fluctuations by a blue surface. A bond-centred 
lattice is introduced and shown in yellow. Long wavelength fluctu¬ 
ations of the order parameter are described by, [Eq. A13]. 


closely follows these references. 

In order to demonstrate the validity of the approach, 
we first derive the continuum theory from a simple spin- 
1 model, ^bbqM^] [Eq- 4]. This is useful, since one can 
check that the resulting action reproduces the results of 
flavour-wave theory [see Section II] at long wavelength. 
Using the results of Section IV, the continuum model can 
be mapped onto the site-centred lattice. Thus we can 
describe the low-energy fluctuations of the spin-nematic 
state found to exist in [Eq. 1]. 

Previous approaches to understanding the spin- 
nematic state at h = 0 in [Eq. 1] reformulate the 

problem in terms of a lattice gauge theory, and solve this 
using a large-N mean-field approach 33 35 . This gives, in 
principle, all the excitations of the spin-nematic state. As 
such, it is interesting to compare the h = 0 field theory 
derived in this Appendix with the lattice gauge theory 
calculations. 


1. Deriving the action 

Before embarking on the calculation, we briefly sum¬ 
marise the mains steps, following the same logic as in 
Refs. [21,23]. First, we note that in [S] [Eq. 4] there 
is an SU(3)-symmetric point at J 22 = 0. When deriving 
a continuum theory, we consider small, but non-zero J 22 , 
and make an expansion around the high-symmetry point. 

The starting point of the calculation is the simplest 
subunit of the lattice, the square plaquette shown in 
Fig. 17. The spin-1 degrees of freedom on this plaquette 
are described using d-vectors [Eq. 7], and the energy is 
given by 77^ [d] [Eq. 9]. The d- vectors obey a length 
constraint d • d = 1 and the phase is set using d 2 = d 2 . 
Energy is minimised on a plaquette by selecting real d 
vectors, with orthogonal alignment between 1st neigh¬ 


bours and parallel alignment between 2nd neighbours. 
A set of matrices can be defined that act on the four d- 
vectors, and allow any configuration to be reached. When 
these act on a ground-state configuration of the plaque¬ 
tte of d- vectors, some of the matrices leave the energy 
invariant, others lead to an energy increase proportional 
to J 22 , while still more lead to an energy increase propor¬ 
tional to Jn, J 12 or a combination of the two. The first 
two sets of matrices will become the low-energy modes in 
the continuum theory, while the third set are considered 
high-energy modes and are eliminated by a Gaussian in¬ 
tegral. 

The square plaquette defines a basic unit from which to 
build a square lattice, and continuum fields are defined at 
the centre of plaquettes. Assuming that the system has 
at least local 2-sublattice AFQ order, the plaquettes can 
be stitched together and an action derived in terms of the 
continuum fields. This action also includes a dynamical 
term, arising from the quantum mechanical overlap of 
nearby director configurations. 

One choice for the ground state of a plaquette is given 

by, 


d| s = (1,0,0), d| s = (0,1,0), (Al) 


where A and B label the two sublattices [see Fig. 17]. In¬ 
cluding fluctuations, and assuming the plaquette is close 
to a ground state configuration, the d vectors can be 
approximated by, 


d A ~ U(0) 

d B « U(0) 
Here, 
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are the relevant subset of SU(3) generators. The param¬ 
eters Zf, l\, v\ and v\ describe small cantings of the pla¬ 
quette configuration away from the low-energy subspace, 
and have been included in Eq. A2 at quadratic order. 
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Plaquettes are stitched together to form a lattice by 
promoting the parameters 0, l \, l \, and r>g to fields, 
and these are defined at the centres of the plaquettes. 
These fields are allowed to fluctuate in time and space, 
and the continuum limit is taken under the assumption 
that the plaquette configurations only vary significantly 
over long lengthscales. 

The action is given by, 

<Snl<Tm = \J [C^ am + £n| CTm ] , (A5) 

where, 

m = 2d A .<9 T d A + 2d B .<9 T d B , (A6) 

and, 

^nlam = (^)plaq? (^7) 


where ( / H) p i a q refers to the Hamiltonian on a single pla¬ 
quette. After substituting in Eq. A2 for the d vectors, 
expanding to quadratic order in the high-energy canting 
fields, making a gradient expansion around the plaquette 
centres and rewriting, 


U(r,r) 


^ n A (r, t) n x B (r,r) n£(r,r)\ 
n\(r,T) n y B (r, r) n y c { r,r) , 

\n z A (r,T) n|(r,r) n z c {r,r) ) 


(A8) 


one arrives at, 

£n'i"m = 2n A d T n A + 2n B <9 T n B + 4^(n A <9 T n B - n A <9 T n B ) 
- 4^2(n A 3 T n B + n A 3 T n B ) + 2w A (n c d T n A + n c <9 T n A ) 

+ 2zt>|(n B d T n c + n B d r n c ), (A9) 


and, 

>£nlcrm = 32Jn ((/i) 2 + (^) 2 ) + 32J 22 (^) 2 
+ I6J22 ((v A ) 2 + (w|) 2 ) - 4J22 (n A n A + n|n|) 

+ E { 4 Jn|n A d A n B | 2 

A=x,y 

+ 2J 22 [( 5 a n A ) 2 + (d A n A ) 2 + (d A n B ) 2 + (d A n B ) 2 ] 

+4Ji2 [| 9 A n A | 2 + |c>a n B | 2 - |n A < 9 A n A | 2 - |n B d A n B | 2 ]} . 

(A10) 


The fields iia and ne inherit the length and phase con¬ 
straints of the d vectors and are further constrained to be 
orthogonal, h/\ ne = 0. The auxiliary field nc = nA x iib 
is introduced as a convenient piece of book-keeping, and 
is not an independent degree of freedom. 

The high-energy, canting fields can be eliminated by a 
Gaussian integral, giving, 

l\ = - 1 (ri A (9riiB - n A <9 r n B ), 

ioJii 

li= 16 (J„ + A,) (a<arnB + n<arfiB) 

'i 

(ncgriiA + n c <9 T n A ) 

lo J 22 

v b = -T^-f—(n B d T nc + n B <9 T n c ). (All) 

loJ 22 


It follows that the action is, 


<Sni CTm [n A , n B , n A , n B ] ~ ^ J dr J d 2 r{ 
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+ 2J 22 [(<9 A n A ) 2 + (S A n A ) 2 + (<9 A n B ) 2 + (<9 A n B ) 2 ] 

+ 4Ji2 [|n A | 2 + |<9 A n B | 2 - |n A <9 A n A | 2 - |n B <9 A n B | 2 ]}. 

(A12) 


This action can be rewritten by reintroducing the matrix 
U(r,r) [Eq. A8] as, 


<Snl< 7 m[U] = ^ J dr J d 2 r { 
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where, 
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This formulation makes clear the SU(2) symmetry of the 
order parameter. In <S n icrm[U] [Eq. A13] we have in¬ 
troduced a set of hydrodynamic parameters. The rela¬ 
tionship between these and the exchange parameters of 
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5 nlcrm [U] [Eq. A13] 

n^[S] [Eq. 4] 

Q,z 

Xo 

[4( Jll + eT22 ) ] 1 

S,z 

Xo 

[4Jn] _1 

Xo Y 

[8 J22] _1 

Po ,Z 

2( Jn + 2Ji2 + 2 J 22 ) 

Po Z 

2( J 11 + 2Ji2 — 2 J 22 ) 

p7 

2( J 11 + J 22 ) 

Aq z 

8V J 11 J 22 


TABLE II: Dictionary for translating between the parameters of 
the continuum field theory for 2-sublattice AFQ order, <S n io- m [U] 
[Eq. A13], and the microscopic model [Eq. 4]. 


^bbqM^] [Eq. 4] is shown in Table II. These can be com¬ 
pared to the hydrodynamic parameters presented in Ta¬ 
ble I. Using the relation v 2 = p/%, one can see that these 
are equivalent at h = 0, except for the velocity (vq’ z ) 2 , 
which differs by a factor proportional to Jf 2 , which is 
considered a small parameter in this Appendix. 


From this one can calculate the dispersion relations of 
the 4 modes as, 
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where Co£ y 0 is twofold degenerate in the 2-sublattice Bril- 
louin zone. This describes 3 Goldstone modes, cj£ y 0 and 
co^q, and 1 gapped mode, cj^’q. 

The Goldstone mode cu® q corresponds to real rotations 
of the d vectors in the ordering plane. This in turn results 
in rotations of the quadrupolar order parameter in the 
ordering plane. 

The pair of Goldstone modes described by co£ y 0 corre¬ 
sponds to real rotations of the d vectors out of the order¬ 
ing plane. This results in rotations of the quadrupolar 
order parameter out of the ordering plane. 

The gapped mode o corresponds to an imaginary 
rotation of the d vectors in the ordering plane. This 
results in spin fluctuations perpendicular to the ordering 
plane, which are out of phase on the 2 sublattices and 
can be thought of as a dynamical spin-density wave. 


2. Linearising the action 


3. Imaginary part of the dynamical spin 
susceptibility 


It is useful to linearise <S n icrm[U] to further bring out 
the physical content of the action. At linear order one 
can approximate, 


U(r, t) 


1 

—01 + i<p 4 

V 02 


01 + 204 — 02 \ 

1 03 I > 

-03 1 ) 


and therefore, 


(A15) 


Inelastic neutron scattering measures the imaginary 
part of the dynamical spin susceptibility, 7s7nx a/3 (q, cj) 
[Eq. 3]. This can be calculated from the linearised ac¬ 
tion, <S n |am[0] [ E q. A16]. 

In order to do this it is first necessary to determine how 
the spin moments are related to the quantum fields. This 
can be achieved by dividing the d- vectors into a real and 
imaginary part, d = u+iv, and noticing that S = 2uxv. 
To lowest order in each of the spin components one can 
show 21 ’ 23 , 


<$n 


lorn 




Xo’ Z (<9t0i) 2 + J2 Po ’ Z ( d *0 if 

A=x,y 


+ Xo’ Z ( < 5t04) 2 + Po’ Z (^A04) 2 + Xq’ Z (^o’ Z ) 2 04 

A=x,y 


+ Xo Y [(^r02) 2 + (9 r 03) 2 ] + Po y [(^A02) 2 + (^A03) 2 ]}- 

A=x,y 

(A16) 


S A = 

^ -20 2 0 4 + 2Xo 4>idt4>2 + Xo’ z 4>2dt(t>i + Xo’ z Xo 

2x7 d tfa 

\ 20 4 - X?’ z <A01 

S B = 

( -2x X odt4>3 

-20304 + 2xo y 0i<9 t 03 - Xo’ Z <t>3d t (pi - Xo’ Z Xo y dtfodtfa 
V 2 04 - Xo’ Z 9t<P l 

(A18) 
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FIG. 18: (Color online). Predictions for the imaginary part of the dynamical spin susceptibility for a 2-sublattice, spin-1 antiferro- 
quadrupolar (AFQ) state in the absence of magnetic field. ^ rn X^\ am an d ^ rn X^ crm [Eq- A20] are calculated from S n \(j m [(f)] [Eq. A16]. 
Hydrodynamic parameters are taken from Table II with Jn = 1, J 12 = 1 and J 22 = 0.5. Dashed red lines show q, CL;* y 0 and All 

predictions have been convoluted with a gaussian to mimic experimental resolution. The circuit T-X-M-T in the bond-centred Brillouin 
zone is shown in Fig. 5. The same linear, normalised colour intensity scale is used as in Fig. 6. 



q q 

FIG. 19: (Color online). Predictions for the imaginary part of the dynamical spin susceptibility for a 2-sublattice, spin-1/2 antifer- 
roquadrupolar (AFQ) state in the absence of magnetic field. The mapping described in Section IV is used to transform and 

^ m XnL-m [Eq- A20] into the site-centred Brillouin zone. An arbitrary set of hydrodynamic parameters have been chosen as (v ® ,z ) 2 = 48, 
(^q Y ) 2 = 12, (r?Q ,z ) 2 = 16 and (Aq’ z ) 2 = 32. Dashed red lines show q, a;* y 0 and q. All predictions have been convoluted with a 
gaussian of FWHM = xxx to mimic experimental resolution. The same linear, normalised colour intensity scale is used as in Fig. 6. The 
results in this Figure can be compared directly with the predictions of the lattice gauge theory given in Ref. 33—35. 


It follows that, 

9TO XH<rm(q- W ) = 9ffl Xlm(q. W ) + u) = 

Sm Xnlam(q» W ) = |xo’^q.O“ W q,o)> ( A19 ) 

and, 

^"ixLm(qM +q,w) = 

Sm xXm(qM + q, W) + + q, w) = 

7r X X o y <V( w -<o) 

^"IXnlam (qM + q, w) = s f 5z % - Wqfo), (A20) 
Xo’ w q ; 0 

where q « 0. This is shown in Fig. 18, and can be 
compared to Fig. 7f, where the physics can be seen to be 
qualitatively the same, despite the different values chosen 
for the parameters Jn, J 12 and J 22 . 

In order to study the bond-nematic phase found in 
[Eq. 1 ], the mapping described in Section IV can 


be used. The field theory thus describes the dynamic 
spin susceptibility close to q = 0, and this is shown in 
Fig. 19. 

4. Comparison with lattice gauge theory 

Finally it is interesting to compare the continuum 
model developed in this Appendix with previous work 

5_l/2 

studying the h = 0 spin-nematic region of _j 
[Eq. 1 ] 33-35 . 

In Ref. [33] a matrix-formed action for the gapless 
modes is written down [Eq. 46 of this reference]. Compar¬ 
ing this to Eq. A13, one can see that these actions have 
the same low-energy form if the generator A 4 [Eq. A4] 
is ignored. The relationship between the hydrodynamic 
parameters of Eq. A13 and Ref. [33] is given by, 

X?’ Z = 2c 2 , Xo y = C 2 , Po’ Z = 2 c 3, 

XV , 

Po = °1 + c 3, C 3 = C4. 


(A 21 ) 
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Comparison can also be made to Ref. [35] in which the 
imaginary part of the dynamical spin susceptibility is cal¬ 
culated via a 1/N expansion scheme. At small q the field 
theory developed in this Appendix agrees well with this 
1/N expansion scheme, which can be see by comparing 
Fig. 19 with Fig. 9 and Fig. 10 of Ref. [35]. In both cases 
there are three Goldstone modes at q = 0, a degenerate 
pair (at low energy) that appear in and a third in 


the longitudinal channel ^mx zz . In Ref. [35], there also 
appear a number of gapped modes close to q = 0 with¬ 
out significant spectral weight, but it is not completely 
clear if one of these is equivalent to the dynamical spin 
density wave that appears in <S n icrm[U] [Eq. A13]. Finally 
Ref. [35] finds a spinon continuum at high energies, and 
this cannot be captured by the low-energy continuum 
theory presented in this Appendix. 
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